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Abstract 

This report presents the results of an extensive micromechanical modeling effort for 
woven metal matrix composites. The model is employed to predict the mechanical response of 
8-hamess (8H) satin carbon/copper (C/Cu) composites. Experimental mechanical results for this 
novel high thermal conductivity material were recently reported by Bednarcyk et al. (1997) along 
with preliminary model results. 

The micromechanics model developed herein is based on an embedded approach. A 
micromechanics model for the local (micro-scale) behavior of the woven composite, the original 
method of cells (Aboudi, 1987), is embedded in a global (macro-scale) micromechanics model 
(the three-dimensional generalized method of cells (GMC-3D) (Aboudi, 1994)). This approach 
allows representation of true repeating unit cells for woven metal matrix composites via GMC- 
3D, and representation of local effects, such as matrix plasticity, yam porosity, and imperfect 
fiber-matrix bonding. In addition, the equations of GMC-3D were reformulated to significantly 
reduce the number of unknown quantities that characterize the deformation fields at the micro- 
level in order to make possible the analysis of actual microstructures of woven composites. The 
resulting micromechanical model (WCGMC) provides an intermediate level of geometric 
representation, versatility, and computational efficiency with respect to previous analytical and 
numerical models for woven composites, but surpasses all previous modeling work by allowing 
the mechanical response of a woven metal matrix composite, with an elastoplastic matrix, to be 
examined for the first time. 

WCGMC is employed to examine the effects of composite microstructure, porosity, 
residual stresses, and imperfect fiber-matrix bonding on the predicted mechanical response of 8H 
satin C/Cu. The previously reported experimental results are summarized, and the model 
predictions are compared to monotonic and cyclic tensile and shear test data. By considering 
appropriate levels of porosity, residual stresses, and imperfect fiber-matrix debonding, 
reasonably good qualitative and quantitative correlation is achieved between model and 
experiment. 
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1. Introduction 


1.1 Woven Composites 

Research on the manufacturing, testing, and modeling of woven and braided composites 
has increased significantly in recent years. The reinforcement phase of these composites consists 
of a woven or braided fabric formed by individual fibers, or by bundles of fibers, called yams. 
One or more layers of the woven or braided fabric are used to reinforce traditional matrix 
materials. It is interesting to note that the concept of woven composites is not a new one. 
Ancient Egyptians used cotton fabrics impregnated with resin to protect fragile mummies. The 
effort to use woven composites for thermal and structural applications, though possibly less 
captivating, is considerably more recent. 

By incorporating a woven reinforcement phase into a composite, rather than utilizing 
unidirectional fibers only, several benefits are realized. A single ply of a woven composite can 
have equivalent thermomechanical properties in several directions. A single ply reinforced by a 
biaxial weave is geometrically similar to a [0°/90°] laminate, while a ply reinforced by a triaxial 
weave can mimic a [0°/±60°] laminate, as illustrated in Figure 1.1. In contrast, a unidirectional 
ply often has poor thermomechanical properties transverse to the fiber direction due to the lack of 
continuous reinforcement in this direction. The deficiency of a continuous ply in the transverse 
direction is often exacerbated by a weak fiber/matrix interface. 



Figure 1.1. a) Biaxial weave pattern; b) Triaxial weave pattern (Chou, et al., 1986). 
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Woven and braided composites usually have superior out-of-plane properties with regard 
to impact and crack resistance relative to composites laminated with unidirectional plies. In fact, 
coated fabrics, which are in essence woven composites, are used to make bullet proof vests. 
Movement of the reinforcement weave can distribute the energy of an impact throughout many 
yams, and, if a crack does form, there are fibers oriented in at least two distinct directions to 
inhibit crack growth. Since the reinforcement phase has a tendency to remain intact 
independently from the matrix, woven composites are less prone to delamination and splitting 
along the fibers. The work of Kaliakin et al. (1996) indicates that these properties give woven 
composites the potential to serve a major role in concrete structure strengthening through the 
application of woven composite plates or jackets directly to the concrete surface. 

Finally, and perhaps most importantly, a woven or braided reinforcement phase offers 
superior stability during composite manufacturing compared to unidirectional fibers. A weave of 
fiber yams is much simpler to handle than the thousands of individual fibers used in a graphite 
fiber reinforced composite, for example. The benefits of the weave’s dimensional stability go 
further. Preforms with complex shapes can be woven or braided from the fiber yams. These 
shapes can then be infiltrated with a metal or epoxy matrix to form a composite in the shape of 
the preform. This procedure is not feasible if individual unidirectional fibers are used. In 
addition, three-dimensional weaves and braids can be produced (see Figure 1.2). A third 
dimension of reinforcement can improve the properties of the composite even further. 




THREE-DIMENSIONAL cylindrical 
CONSTRUCTION 


THREE-DIMENSIONAL BRAIDING 


Figure 1.2. Examples of 3-D weave patterns, a) Cylindrical construction, 
b) 3-D braiding (Chou, et al., 1986). 


1.2 Modeling of Woven Composites 

While woven composites offer advantages over unidirectional composites and laminates, 
they are also more challenging to model. The woven reinforcement phase consists of yams that 
undulate in and out of a plane. Thus the geometry of the composite is inherently three- 
dimensional. In addition, there are many different ways in which the reinforcement fabric can be 
woven (see Figure 1.3), and each of these weaves has a different repeating unit cell. These 
factors combine to make modeling of woven composites challenging. 



Plain Weave 



2/1 Warp, Twill 




2/2 Cord 



Figure 1.3. Examples of biaxial weaves (Miller, 1968). 


1 .2. 1 Finite-Element and Boundary-Element Models 

The complex three-dimensional geometry of woven composites makes finite-element 
modeling difficult. One can imagine constructing a three-dimensional mesh for finite-element 
implementation for each type of weave. If the geometry is then slightly changed, for example if 
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adjacent yams are placed closer together, an entirely new mesh would be necessary. The effort 
required for such numerical modeling may be prohibitive. However, this type of effort was 
undertaken by Dasgupta and Bhandarkar (1994) and Dasgupta et al. (1996) to model the elastic 
behavior of a plain weave glass/epoxy composite. In this investigation, reasonable elastic 
constants were predicted for a realistic geometric representation of the composite, but only with 
great computational effort. 

Whitcomb et al. (1992) proposed a finite-element model for woven composites in which 
spatial variations of material properties were accounted for within a single element. This 
approach could potentially decrease the number of elements required to accurately model the 
geometry of a woven composite. However, a traditional three-dimensional finite-element 
analysis was used by Whitcomb and Srirengan (1996) to model progressive failure in plain 
weave graphite/epoxy composites with varying degrees of fiber waviness. The finite-element 
approach was also used by Glaessgen et al. (1996) to examine the internal displacement and 
strain energy density fields in a plain weave glass/epoxy composite. Here, geometric 
complexities inherent to woven reinforcements (which greatly affect internal fields) are 
accounted for, but at a high computational cost. 

Marrey and Sankar (1997) performed elastic finite element analyses of plain weave and 
5-hamess satin composite plates through the use of homogenized brick elements. Effective 
properties of the brick elements, which represent the composite repeating unit cell, were first 
determined via finite element analysis. Then these homogenized elements were assembled, 
under appropriate boundary conditions to form a plate. It should be noted that the above finite- 
element analyses considered only woven composites with elastic phases. Due to the complex 
geometry of woven composites (and thus the large number of degrees of freedom), inclusion of 
matrix inelasticity in finite element models for these materials would require immense execution 
times. Analysis of a simple woven metal matrix composite (with an elastoplastic matrix) via 
finite elements may be possible through the use of supercomputers, but to date, such an analysis 
has not been reported. 

A boundary-element model developed by Goldberg and Hopkins (1995) that has been 
used to examine the elastic response of woven composites deserves reference. The boundary- 
element method requires less computational and mesh generating effort than the finite-element 



method, yet it can offer similar geometrical accuracy for woven composites. A version of this 
model with matrix inelasticity is under development and may have potential for modeling woven 
metal matrix composites (Goldberg, 1996). 


1 .2.2 Approximate Analytical Models 

Another route to modeling woven composites was taken by Chou and Ishikawa (1989). 
These authors have developed a well-known series of models based on classical lamination 
theory for predicting the thermoelastic response of certain types of woven composites. The 
mosaic model treats the weave as an assemblage of cross-ply laminates; however, only a two- 
dimensional portion of the actual repeating unit cell of the composite is considered (see Figure 
1 .4). As indicated, the cross-ply sections are assembled under an iso-stress or iso-train condition. 
Similarly, the crimp model considers only a two-dimensional portion of the actual repeating unit 
cell, but adds crossover of the fiber yams in the loading direction (see Figure 1 .5). The bridging 
model combines the mosaic and crimp models (see Figure 1 .6) by taking weighted averages of 
effective stiffness terms in an attempt to account for the three-dimensionality of the actual 
repeating unit cell. This model still represents a highly idealized geometric representation of a 
woven composite, but, like the mosaic and crimp models, it offers the ability to model 
composites reinforced with some more complex weaves beyond the plain weave pattern. 

Naik and co-workers have developed several models for plain weave composites based 
on the approach of Kabelka (1980, 1984). Kabelka's approach (which was also originally applied 



Figure 1.4. Mosaic model geometry. 
(Chou and Ishikawa, 1989) 



Figure 1.5. Crimp model geometry. 
(Chou and Ishikawa, 1989) 
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Figure 1.6. Bridging model geometry (Chou and Ishikawa, 1989). 

to a plain weave composite) considers a repeating element from a single cross-section of the 
composite in each of the yam directions. The geometry is shown in Figure 1.7, where warp and 
weft refer to the two fiber directions in the biaxial plain weave reinforcement pattern (weft is 
commonly referred to as fill). Expressions were developed for local dimensions and fiber angles 
in these cross-sections. The effective properties of the cross-sections were then taken as the 
mean integral values of the local effective properties determined via classical lamination theory. 
This method is unrealistic in that it models the behavior of the entire woven composite as the 
behavior of one cross-section or “slice” of the geometry in each direction, while in reality the 
cross-sectional geometry is changing throughout the composite. Naik and Ganesh (1992) 
remedied this limitation by performing a Kabelka-type analysis on a number of “slices” from the 
actual three-dimensional plain weave composite unit cell, and assembling these slices under an 
iso-strain condition. The authors refer to this model as the Slice Array Model (SAM) (see Figure 
1.8). The Element Array Model (EAM) considers slices taken in both in-plane directions such 
that discrete elements, rather than slices, are formed (see Figure 1 .9). Each element is modeled 



WARP 


Figure 1 .7. Cross-section geometry modeled by Kabelka (1984) 



UNIT CELL 




Figure 1.9. EAM model geometry (Naik and Ganesh, 1992) 
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with classical lamination theory, and the elements are assembled in one in-plane direction to 
form slices, and then in the transverse in-plane direction to form the actual repeating unit cell of 
the plain weave composite. The elements are assembled under the iso-stress condition along the 
loading direction and under the iso-strain condition transverse to the loading direction. The order 
in which these two assembly processes proceed distinguishes two distinct models whose 
predictions can vary significantly. Comparison with experimental in-plane elastic modulus data 
for plain weave graphite/epoxy shows that one or both models are reasonably accurate for the 
various composite properties. A similar iso-strain approach was employed by Naik (1995) to 
develop a general code for elastic analysis of woven and braided composites called TEX-CAD. 
This analytical model allows analysis of a wide range of geometries and includes damage 
accumulation, composite failure, and yam bending. 

Another analytical model, developed by Karayaka and Kurath (1994), uses a 
homogenization technique in conjunction with classical lamination theory. Effective 
(homogeneous) properties of a single ply of a woven composite representative volume element 
are determined via a unit cell analysis in which all in-plane strain components and out-of-plane 
stress components are assumed to be constant throughout the composite. The effective properties 
of the woven composite plies are then used in classical lamination theory to model a nine-ply 5- 
hamess satin weave graphite/epoxy laminate. 

Thus, it is clear that a considerable amount of effort has been expended in an attempt to 
model woven composites analytically. Most models have been shown to be reasonably 
successful at predicting the effective elastic properties of woven polymer matrix composites. 
However, like numerical models, all analytical models for woven composites reported to date 
lack the ability to simulate the inelastic constitutive behavior of metal matrix composites. 

1.3 Objectives of Present Investigation 

The primary objective of this investigation was to develop an analytical model for woven 
metal matrix composites which is both realistic and practical. Previous models for woven 
composites do not incorporate inelastic behavior of the matrix and thus are insufficient for woven 
metal matrix composites. The present model is based on an embedded approach in which a local 
model is embedded in a global model. The global model is an extension of Aboudi’s (1994) 
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three-dimensional generalized method of cells (GMC-3D). This model simulates the overall 
behavior of the woven composite through analysis of the actual three-dimensional repeating unit 
cell. The local model is an extension of Aboudi’s (1987) original method of cells. This model 
simulates the microscopic behavior of the woven composite, on the level of the individual fibers 
and matrix which constitute the infiltrated fiber yams in the woven composite. 

The extensions of both models mentioned above are necessitated by the two requirements 
placed on the analytical model, namely that it be realistic as well as practical. From a realistic 
stand-point, a finite element model for 8H satin C/Cu would probably be most desirable. 
However, the large number of three-dimensional elements that this approach would require, 
coupled with cyclic loading, temperature dependence, and matrix inelasticity make finite element 
analysis quite impractical. Conversely, approximate analytical models previously developed for 
woven composites are quite practical, yet they leave much to be desired from the stand-point of 
realism. The analytical model developed herein splits the difference between these approaches. 
It is considerably more practical than finite element models (yet somewhat less realistic) and 
considerably more realistic than previous analytical models (yet somewhat less practical). 

To achieve the above claims of both practicality and realism, it was necessary to modify 
and extend both GMC-3D and the original method of cells. In its original form, GMC-3D is not 
computationally efficient. Preliminary modeling work for woven composites (Bednarcyk et al., 
1997) indicated that this issue becomes problematic as the number of subcells analyzed by the 
model becomes large (as it does in the case of woven composites). To remedy this, the equations 
of GMC-3D were completely reformulated, taking advantage of the constancy of certain stress 
components in certain subcells. This reformulation follows that performed for the two- 
dimensional version of the generalized method of cells (Pindera and Bednarcyk, 1997), and it 
improves the computational efficiency of the model dramatically. In addition, in order to 
accommodate the incorporation of the local model, GMC-3D was extended to include fully 
anisotropic elastoplastic constituents which result from rotating local quantities to the global 
coordinates. 

The original method of cells is embedded within GMC-3D to model the local behavior of 
each three-dimensional subcell. In order to make the analytical model as realistic as possible, the 
original method of cells was extended to include matrix plasticity, imperfect fiber-matrix 
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bonding, and consistent (Brayshaw) averaging to allow the model to simulate the transversely 
isotropic behavior of a unidirectional composite (Brayshaw, 1994). While each of these features 
had been previously incorporated in the original method of cells independently, they had not 
previously been combined. 

Thus, an analytical model for woven metal matrix composites is developed which is both 
realistic and practical. This model is called Woven Composites Generalized Method of Cells 
(WCGMC). The model is versatile and can simulate a wide range of thermomechanical loading 
conditions for woven and braided composites. Herein, WCGMC is employed to simulate the 
mechanical behavior of 8-hamess (8H) satin carbon/copper (C/Cu). This woven metal matrix 
composite is a candidate for high heat flux applications. An extensive experimental investigation 
into the mechanical behavior of 8H satin C/Cu was performed by Bednarcyk et al. (1997), the 
results of which are summarized in Chapter 4 of this report. A parametric study is performed 
with the model to highlight the effects of geometric unit cell refinement, fiber volume fraction, 
porosity, residual stresses, and imperfect fiber-matrix bonding on the mechanical response of 8H 
satin C/Cu. Model predictions are then compared with the experimental data for the composite 
to evaluate the accuracy of the model and identify areas for improvement. 
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2. Analytical Model - WCGMC 

The analytical micromechanics model developed in this investigation is called woven 
composites generalized method of cells (WCGMC). It is based on an embedded approach in 
which a local micromechanics model, the original method of cells, is embedded in a global 
micromechanics model, the three-dimensional generalized method of cells (GMC-3D). This is 
shown schematically in Figure 2.1. GMC-3D uses an arbitrary number of homogeneous three- 
dimensional subcells to represent the three-dimensional repeating unit cell of a material. In 
WCGMC, the equations of GMC-3D are used, but the three-dimensional subcells are permitted 
to be heterogeneous and thus able to represent a portion of an infiltrated fiber yam. This subcell 
heterogeneity is accomplished by using the original method of cells to model the local behavior 
of the three-dimensional subcells. In addition, matrix plasticity, Brayshaw averaging (to allow 
the infiltrated yam subcells to exhibit transversely isotropic behavior) (Brayshaw, 1994), and 
imperfect fiber-matrix bonding are incorporated in WCGMC on the local level in the embedded 
original method of cells. In order to provide WCGMC with the level of computational efficiency 



Figure 2.1. Schematic representation of the incorporation of composite subcells into GMC-3D 

via the original method of cells. 
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necessary to simulate the thermo-inelastic behavior of a sufficiently refined woven composite 
unit cell geometry, it was necessary to reformulate the basic equations of GMC-3D. The 
resulting global equations used in WCGMC represent a reduction of nearly sixteen times in the 
number of unknown global quantities that must be determined in the model. The original 
formulation of the GMC-3D equations will be presented first, followed by the reformulation, and 
finally by the local equations for the subcell behavior. 

2.1 GMC-3D - Original Formulation 

For a complete derivation of the equations for GMC-3D, see Aboudi (1994). The 
geometry considered by the model is shown in Figure 2.2. A multi-phase material is represented 
by a parallelepiped unit cell which repeats infinitely in the three mutually orthogonal directions. 
The cell is divided into an arbitrary number of parallelepiped subcells, each denoted by the three 
indices (a py) • The total number of subcells in each direction is denoted by N a , Np, and N y . 



Figure 2.2. GMC-3D geometry. 
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The displacement field in each subcell is assumed to be linear in the local subcell 
coordinates, (x, , x 2 ^ xj^), centered in the middle of each subcell, 

u WPr) = Jflfir) + + x[ r) if/^ Pr) , (2.1) 

where the subcell microvariables, ^ determine the displacement field 

dependence on each subcell coordinate. The subcell strain components are given by the usual 
kinematic relations, 


4 ^ 4 (^ 45 *’) 


i,j~ 1 , 2 , 3 . 


( 2 . 2 ) 


Because the displacement field is assumed to be linear, equation (2.1), the strain components 
within each subcell are constant, and the average strain components for the unit cell are given by 
the volume average of the strain components in the subcells, 


1 


£,J dhl 


d h l e {apr) 

u a n p l y b i] 


(2-3) 

a P Y 

Note that in order to simplify notation, summations over the indices (a/3y) will be expressed as 
in equation (2.3) above. That is, 


N a Nfi 

Z«Z*. z-z*. I- 1' 

P p= 1 y 


(2.4) 


0=1 p 0=1 r 

Since the subcell strain components are constant within each subcell, the stress 
components in each subcell are constant as well. The subcell stress components are related to the 
subcell strain components by the subcell constitutive equations, 


T (<*Pr) _ (+»)(#>) - £ 


p(aPr) 

Id 


-a ( f r) AT 


)■ 


(2.5) 


where are the subcell plastic strain components, a { ^ Pr) are the subcell coefficients of 

thermal expansion (CTEs), and AT is the temperature change from a reference temperature. The 
average stress components in the unit cell are given by the volume average of the subcell stress 
components, 


*.-£zzz«w-r>. 


( 2 . 6 ) 
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Continuity of displacements and tractions is required between subcells within the unit 
cell, and between the unit cell and adjacent unit cells. These continuity requirements are 
imposed in an average sense, that is the integrals of the appropriate displacement and traction 
components along the appropriate boundaries are required to be continuous. Imposing the 
displacement continuity conditions (see Aboudi (1994) for details) gives rise to six continuum 
equations of the form, 


a 

P= Iv- 


(2.7) 

Yh M?) =hs 
P 

or = 1,.. 

. ,A^ a , y \,... t N y t 

(2.8) 

V/ = 77 

3 * 6 33 > 
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a - 1,... 

> 

II 

** 

(2.9) 

II V,4f ’ ■ 
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= hi S 23 , 


(2.10) 

II dJ/S*' 

a y 

= dls n . 

II 

> 

(2.11) 

Il^vif’ 

= dhe a , 

r = l...,Ny. 

(2.12) 


a 0 

These continuum equations form a system of equations which can be expressed as, 

A c e, = Je, (2.13) 

where, e = {£ n £ 22 £33 2 i 23 2 i ]3 2i 12 } , and z s is the 6N a NpN y - order subcell strain 

vector given by e s = {e^ 111 ^ ... e^° N t Nr )} t where each vector consists of the six subcell 
strain components. a g contains cell geometric dimensions only, and it is an 
[N a (Np + N y + l) + Np(N y + 1) + N r ] x 6N a NpNy matrix. J contains cell dimensions, and it is an 

[N a (Np + Ny+l) + Np(Ny + l)+Ny] x 6 matrix. 

Imposing the traction continuity and using (2.5) gives rise to the system of equations, 

A M (s s -z p s -a s AT)=0, (2.14) 

where z p s and a T s are 6 N a NpN y - order subcell plastic strain and subcell CTE vectors, similar 
in composition to z s . The matrix A w is 6N a N p N y - (N a N p + N a N r + N fi N y 
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(iV a + N p + N r j x 6 N a NpN y and contains the subcell stiffness components. Combining (2.14) 
and (2.13) yields, 


Ae^-D^ + a s ATj = Ke , 


(2.15) 


where, A = 



Equation (2.15) can be solved for e s , 


£ s = A £ + D^8 £ + a S AT ^ , 


(2.16) 


where, A = A 1 K and D = A D . If the matrices A and D are partitioned into N a N pN y sixth- 


order square submatrices such that A = 


aO") 

\Wp n t) 


and D = 


d(' u ) 


, then equation (2.16) 


implies that, 

E w = A (#)g +D (#)j £ p + ttjA r) . (2.17) 

This equation gives the strain components in each subcell in terms of the applied cell strains, the 

subcell plastic and thermal strains, and two concentration matrices, aS°^ and . 
Substituting (2.17) into (2.5) yields, 

*(**) =d^)[ A (^)e + D (^ r) (ef +a s AT) -(s*^ +a (<,fr) Ar)], (2.18) 

and using the average stress equations (2.6), the elements of the effective elastoplastic thermo- 
mechanical constitutive equation, 

a = C*(s -a* AT) , (2.19) 

can be found. The effective elastic stiffness matrix is. 



Z£2Xv, 

a P y 


Q(°Pr) pJ,aPr) 


the cell plastic strain vector is. 


( 2 . 20 ) 
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z p _ e JK^)) t 


( 2 . 21 ) 


the average CTE vector is. 
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a' = -^IISXV, C^D'^a, -«<«»>) , (2.22) 

a is the average stress vector, and e is the imposed average strain vector. 

In this formulation of GMC-3D, the subcell strain components, given by equation (2.17), 
serve as the basic unknown quantities. Since there are six unknown strain components per 
subcell, the total number of unknown quantities which must be determined is 6N a N p N r . 

Determining these quantities (i.e., employing equation (2.17)) involves inverting the 
6N a N fi Nx6N a N fi N y A matrix. In the presence of plasticity, the thermomechanical 

loading must be applied incrementally, and an iterative solution procedure must be employed at 
each load level (see Section 2.6). Thus, for a given thermomechanical loading simulation, the 

unknown subcell strain components must be determined, and the A matrix must be inverted, 
hundreds or thousands of times. As the number of subcells in the repeating unit cell becomes 
large, the original formulation of GMC-3D becomes increasingly computationally inefficient. 
For this reason, the GMC-3D equations have been reformulated so that sufficiently refined 
woven composite unit cells can be modeled. 


2.2 GMC-3D - Reformulation 


Since the individual subcells in WCGMC are heterogeneous, the subcells in GMC-3D 


must be anisotropic. The subcell anisotropic constitutive equations can be expressed as. 
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(2.23) 


The subcell traction continuity conditions require that, at subcell and unit cell interfaces, the 


tractions be continuous. Since the unit cell and the subcells are parallelepipeds, each interface is 


normal to one coordinate axis. Thus, each unit normal vector for each interface is parallel to one 
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coordinate axis, and particular subcell stress components (all of which are constant within a 
subcell) are equal to the traction components at the interfaces. That is, 


(n ) 

t; = VijKj = 


7 n=\ 

a J 

0 rij = 0 


(2.24) 


This allows each traction continuity condition to be expressed in terms of one subcell stress 
component. In fact, the traction continuity conditions which are applicable to the normal subcell 
stress components require that each normal stress component be constant through all subcells 
which are adjacent in the coordinate direction of that subcell stress component. That is, for 

example, is constant when following a row of subcells through the unit cell shown in 

Figure 2.2 in the x ] -direction. This condition can be expressed as, 

J)Pr) _ JlPr) _ _ J s afir) _ Jfir) q 25 ) 

a U -cr ll _y ll > 

where has been introduced to denote the 11 stress component in each row of subcells 
which are adjacent in the x ] -direction. Similarly, for the remaining normal subcell stress 
components, 

°22 r) = °22 2r) = • • • = Cr ( 2 f #y) = 4" r) > (2.26) 

cr<f '> = cr ( 3 f 2) = . . . = a[f Nr } = T^f ] . (2.27) 

The traction continuity conditions which affect the shear stress components can similarly 
be applied. One difference is that, by nature of the symmetry of the stress tensor (i.e. cr . = a jt ), 
two traction continuity conditions affect each subcell shear stress component. For example, 
cr is constant when following a row of subcells which are adjacent in the x 2 -direction, while 

cr 32 ^ is constant when following a row of subcells which are adjacent in the x 3 -directions (see 

Figure 2.2). However, since cr\°f y) = (J'* 2 Pr \ the 23 subcell stress component must be constant in 
each layer of subcells with a constant value of a . This condition can be expressed as, 


r (° 1 ?') _ /T ( a2 r) _ _ A°^f r ) 


0 _ _(« 02 ) _ _ Jpfr) \ 23 


= Cr 32 = 


Safir) _ rA a fr) - r( a ) 

*23 — ° 32 “ y 23 » 


(2.28) 


where has been introduced to denote the 23 stress component in each layer of subcells 
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which are adjacent in the x 2 -direction or the x 2 -direction. Similarly, for the remaining subcell 
shear stress components, 
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(2.29) 


(2.30) 
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The utility of accounting for traction continuity in this explicit manner is clear. There are 
only N fi N r +N a N y + N a Np + N a + Np + N r unique subcell stress components, which have 


been denoted Thus, if these subcell stress components, rather than the 6N a N J3 N y subcell 

strain components, are employed as the basic unknowns, the number of unknown quantities is 
reduced substantially. This reduction in unknowns results in greater computational efficiency for 
the model. 

Substituting for the subcell stress components in the subcell constitutive equations (2.23) 
using equations (2.25) - (2.30), solving for the subcell total strains, and substituting into the 
continuum equations (2.7) - (2.12) yields. 


1 ''ll + 2.^12 '22 + '33 ' + 

a a a 
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(2.32) 


Y i cMi')r(fr) , y / c( a ^) r M Y/ cWrW+V/ . Y/ 
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(2.33) 


a = fi = 1,...,A^ 
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r = l...,N y 


(a0r) 

(2.36) 


These equations can be assembled into a globed equation in matrix form and written, 

GT = f m -f t AT -f p , (2.37) 

where G is an NpN y + N a N r + N a Np + N a + Np + N y -order square matrix containing subcell 
dimensions and subcell compliance components, T is an N p N y + N a N r + N a N p + N a + N p + N y - 

order vector containing the unknown subcell stresses, f" 1 is an 
NpNy + N a N y + N a Np + N a + Np + N y -order vector containing cell dimensions and global (cell) 

strains, f * is an NpN y + N a N y + N a Np + N a + Np + N y - order vector containing subcell dimensions 
and subcell coefficients of thermal expansion, and f p is an N p N y + N a N y + N a Np + N a + Np + N y - 
order vector containing subcell dimensions and subcell plastic strain components. The structure 
of the G matrix is shown in Figure 2.3. It consists of 36 submatrices, only 12 of which are fully 
populated. The components of G , T , f m , f * , and f p are given in the appendix. 
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Figure 2.3. Structure of the global matrix G . The case shown is for N a = 4,N„=2,N r =3 
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To obtain explicit expressions for the subcell stress components, the global equation 
(2.37) is inverted to obtain the subcell mixed concentration equation, 
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, (2-38) 


where the A(*\ B(*\ X^, A ( * } , Q^* } , , rj*\ and 0(* } terms are given in the appendix. 

This equation is referred to as a mixed concentration equation because the local or subcell 
stresses are related to the global or cell strains. The average stress equation (2.6), can now be 
applied to the subcell stresses in equation (2.38), yielding. 
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The expressions for 7^ from equation (2.38) are substituted into equations (2.39), and the 
results compared with the global (cell) constitutive equation, 
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to yield closed-form expressions for the cell effective stiffness components, C* , the cell effective 
coefficients of thermal expansion, a’, and the cell effective plastic strain components, ££ . 
These expressions are, 
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(2.43) 


Thus, as was the case in equations (2.20) - (2.22) in the original formulation of GMC-3D, 
equations (2.41) - (2.43) provide closed-form expressions for the effective thermoelastic 
constants and effective plastic strain components for the three-dimensional unit cell. However, 
comparing equations (2.17) and (2.38) shows that, for applied thermomechanical loading, there 
are far fewer unknown variables to be determined using the reformulated version of the GMC-3D 
equations. This is clearly illustrated in Figure 2.4. The number of unknown variables is plotted 
versus the number of subcells in the repeating unit cell to be analyzed for the case in which the 
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number of subcells is identical in each direction. For 1,000,000 subcells, there are nearly 200 
times fewer unknowns in the reformulated version compared to the original version. 

2.3 Heterogeneous Subcells via the Reformulated Original Method of Cells with 
Imperfect Fiber-Matrix Bonding 

In order to model a woven composite with GMC-3D, it is necessary for the three- 
dimensional subcells to represent the infiltrated fiber yams. This is achieved by allowing the 
three-dimensional subcells to be heterogeneous, with unidirectional fibers. In WCGMC, the 
fiber direction of the unidirectional composite comprising each three-dimensional subcell is 
arbitrary, as are the fiber material, matrix material, fiber-matrix debonding parameters, and fiber 
volume fraction. The unidirectional composite in each subcell is modeled in its principal 
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material coordinates by the original method of cells, also developed by Aboudi (1987) (see 
Figure 2.1). 

One feature of the original method of cells is that it simulates a unidirectional composite 
consisting of an isotropic matrix and a transversely isotropic fiber as cubic rather than 
transversely isotropic. That is, the effective stiffness matrix of such a composite as determined 
by the original method of cells has six independent components when it should have five. In 
order to model the unidirectional composite as transversely isotropic, Aboudi (1987) suggested 
averaging the effective stiffness components in the transverse plane of symmetry. This 
procedure indeed provides transversely isotropic elastic behavior, but Brayshaw (1994) showed 
that it leads to an inconsistency in the unit cell stress components. That is, the weighted sum of 
the subcell stress components is not equal to the average stress components in the unit cell. The 
correction of this inconsistency is imperative when the original method of cells is embedded in a 
global micromechanics model because the average stress components for the original method of 
cells unit cell are continuously exchanged between the two models. Hence, Brayshaw’s 
averaging procedure, which eliminates the inconsistency was employed in WCGMC and will be 
described in Section 2.4. 

Since C/Cu composites are known to exhibit imperfect fiber-matrix bonding, this feature 
was included in WCGMC on the local level in the original method of cells. The original method 
of cells equations for composites with imperfect bonding were provided by Aboudi (1988), and 
the consistent formulation for the original method of cells equations with perfect bonding was 
provided by Brayshaw (1994). However, the combination of imperfect bonding and the 
consistent formulation using the equations provided by the aforementioned authors proved 
unmanageable. Thus, a different approach was taken. The equations of the original method of 
cells were reformulated along the lines of the GMC-3D reformulation. This original method of 
cells reformulation amounts to a specialization of the reformulation of GMC-2D presented by 
Pindera and Bednarcyk (1997), with the additional features of imperfect bonding and Brayshaw 
averaging. As will be shown, the reformulated equations of the original method of cells lend 
themselves to the inclusion of imperfect bonding and consistent averaging. 

The original method of cells provides effective constitutive equations for the two- 
dimensional rectangular unit cell shown in Figure 2.5 consisting of three rectangular matrix 
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Figure 2.5. Original method of cells geometry. 

subcells and one rectangular fiber subcell. Each subcell in the original method of cells is denoted 
by the two indices ( ffy ). These constitutive equations, in turn, describe the average response of 
the subcells in GMC-3D. It is important to distinguish between the subcell indicial notation used 
in the original method of cells (py) and that used in GMC-3D (apy ). Though the notation is 
similar, the indices refer to the subcell quantities in two distinct models. The procedure for 
developing the effective constitutive relations in the original method of cells is similar to the 
procedure used in GMC-3D: a first order displacement field is assumed for the subcells and 
continuity of displacements and tractions between subcells and between cells is imposed in an 
average sense. 

Imperfect fiber-matrix bonding in the original method of cells is modeled by introducing 
a discontinuity in the displacement components at the fiber-matrix interface. This displacement 
discontinuity allows slippage or separation at the interface, and it can be characterized by a linear 
relationship with the traction components at the interface (Aboudi, 1988). The discontinuity in 
the displacement component normal to a given interface is related to the normal traction 
component at the interface, while the tangential displacement discontinuity is related to the 
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tangential traction component at the interface. These two relations can be written, 

[«„],- *,<%!,. C- 44 ) 

[«,!, = (2.45) 

where [w n ] ; and [u,^ are the discontinuities, or “jumps” in the normal and tangential 

displacements at the interface /, and R n and R, are the normal and tangential debonding 
parameters. These parameters are, in effect, compliances of a flexible fiber-matrix interface, and 
R n -*■ 0 , R, 0 corresponds to a perfectly bonded interface, while R„ -> oo , /?, — > co 
corresponds to a completely debonded interface. 

In the original method of cells, only subcell interfaces between the fiber and matrix are 
permitted to have imperfect bonding. The debonding parameter notation which will be used is: 

• => applies to interface between subcell (fi,y) and subcell (/? + \,y ) , 

• R^ ri => applies to interface between subcell (/?,y) and subcell (/?,/ + 1) . 

In R-lf 3 ^ , i and j determine which debonding parameter is referred to, thus. 


={*” or ! ? i=J >. 2 * 1.1 . 

J 1 7?, or 0 for / * j 


(2.46) 


Note that an debonding parameter is disallowed since the displacement in the x, -direction 
is required to be uniform throughout the unit cell. 

In this section, for notational simplicity, summation over the indices f3 and y will imply 
a range of 1 to 2, since the repeating unit cell consists of two subcells in each direction (see 
Figure 2.5). That is, 

2> = i> !• = £•. <2.47) 

P P= 1 y r« 1 

As was the case for the GMC-3D equations, application of displacement continuity 
results in a set of continuum relations (see Aboudi (1988) and Sankurathri et al. (1996) for 
details). These continuum relations now contain the debonding parameters and can be written. 


+ R 22 ] )<J 22 ') = h£ 22 . Y = h 2 , ( 2 ‘ 48 ) 
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The normal constitutive equations for the homogenous subcells are, 
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where, e n is the imiform axial strain for all subcells, and s[f^ are the remaining normal 
subcell strain components, are the subcell CTEs, are the subcell plastic strain 

components, Sjf r ' > are the subcell compliance components, and are the 

subcell stress components. Solving this system of equations (2.50) for af^ yields, 

jfir) — _UV. .jfir), (2 5n 

CT 11 " a H *11 **12 ct 22 **13 °33 )■ 


Substituting this expression (2.51) in the second and third subcell constitutive equations (2.50) 
yields. 
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Substituting these equations (2.52) and (2.53) into the two continuum equations (2.48) and (2.49) 
yields, 
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As was the case in GMC-3D, imposing continuity of normal tractions amounts to 
requiring normal stress components to be constant in rows of subcells which are adjacent in the 
appropriate normal direction. For the original method of cells, this is simpler since there are only 
four subcells in the repeating unit cell. The normal traction continuity conditions can be 
summarized as. 
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Through the use of these equations (2.56) and (2.57), along with the grouping of terms, equations 
(2.54) and (2.55) can be written as, 
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Examining the terms defined in equation (2.60) and (2.62) (which contain debonding 
parameters), with the use of equation (2.46), it is found that, 
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In matrix form, equations (2.58) and (2.59) can be written. 
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If the components of the inverse of the matrix on the left hand side of equations (2.69) are 
denoted by , these equations can be solved for the subcell stresses, 
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Equations (2.70) and (2.71) relate the subcell normal stresses to the global (applied) strains and 
terms involving plastic and thermal strains. Equations (2.70) and (2.71) are thus normal mixed 
concentration equations for the original method of cells. Substituting equations (2.70) and (2.71) 

into the expression for <?\ Pr K equation (2.51), yields, 
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(2.74) 


where. 
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The stress averaging equations require that the geometrically-weighted subcell stresses sum to 
equal the global (cell) stresses, 

±YLVr°u ) (2.76) 

P Y 

iXI (2.77) 

P Y Y 

41^=533. (2.78) 

/5 r "A 

Substituting the expressions for the subcell stresses, equations (2.70), (2.71), and (2.74), into the 
averaging equations, (2.76), (2.77), and (2.78), and comparing with the form of the global (cell) 
normal constitutive equation. 
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allows the identification of the effective normal stiffness components, C* , coefficients of thermal 
expansion, a* , and plastic strain components, H , 
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Equations (2.79) - (2.82) provide the effective normal constitutive equation, stiffness 
components, CTE components, and plastic strain components needed to model the heterogeneous 
subcells in WCGMC. 

To continue, the shear constitutive relations for the reformulated original method of cells 
will be developed. The shear continuum equations are (see Aboudi (1988) and Sankurathri et al. 
(1996) for details), 
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Substituting for the subcell shear strains in the above equations (2.83) - (2.85) using the 
orthotropic subcell shear constitutive equations. 
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and making use of the subcell shear traction continuity conditions, which can be written, 
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Examining the debonding parameters in the above expressions (2.93) and using equation (2.46), 
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) can be easily solved for the subcell shear stresses, 
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Equations (2.98) - (2.100) relate the subcell shear stress components to the average (cell) strain 
components and plastic terms. Thus, these equations are the shear mixed concentration 
equations for the original method of cells. 

The shear stress averaging equations, 
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can then be applied using equations (2.98) - (2.100) to yield, 
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Comparing these equations (2.107) - (2.109) to the global (cell) orthotropic shear constitutive 


equations, 
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allows the identification of the cell effective shear stiffness components and plastic strain 


components. 
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Thus, closed-form constitutive equations, (2.79) and (2.1 10), have been developed for the 
original method of cells with imperfect fiber-matrix bonding. Note that the debonding 
parameters appear only in the terms defined in equations (2.67), (2.68), and (2.97). Hence, 
unlike the original formulation of the original method of cells, including imperfect bonding in the 
reformulation amounts to simple modifications of a few equations. 

In WCGMC, given the geometry and properties of the fibers and matrix that constitute 
the unidirectional composites in each three-dimensional subcell, the reformulated original 
method of cells is used to calculate effective thermoelastic properties for each three-dimensional 
subcell in the principal material coordinates of the given three-dimensional subcell. An effective 
stiffness matrix is then assembled for each three-dimensional subcell in the principal coordinates. 
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and this stiffness matrix is rotated in three dimensions to the coordinates of the three-dimensional 
cell. The reformulated version of GMC-3D uses these effective stiffness matrices for each three- 
dimensional subcell to determine the effective stiffness matrix and effective properties for the 3D 
cell. If the woven composite is treated as purely elastic, this is all that is required. However, if 
matrix plasticity is included, the formulation becomes more complex. 

The subcell plastic strains are determined from the original method of cells in conjunction 
with the classical incremental plasticity theory (see Section 2.5). The strains in each three- 
dimensional subcell, known from the solution of the global constitutive equation (2.40), the 
mixed concentration equations (2.38), and the subcell constitutive equations (2.23), represent the 
average unit cell strains in the original method of cells. From these average unit cell strains, 
stress and strain components for each of the four subcells can be determined via the original 
method of cells mixed concentration equations and subcell constitutive equations. The 
knowledge of these subcell stresses and strains allows the determination of plastic strains within 
each original subcell using the classical incremental plasticity theory equations presented in 
Section 2.5. These plastic strains are then used to determine global (cell) plastic strains in the 
original method of cells via equations (2.82) and (2.112). These cell plastic strains are then 
rotated back to the global coordinates, and then they represent the three-dimensional subcell 
plastic strains for the three-dimensional subcell considered. The three-dimensional subcell 
plastic strains are used to re-solve the global constitutive equation (2.40), for the entire three- 
dimensional unit cell. 

2.4 Brayshaw Averaging 

When the original method of cells is used to model a unidirectional composite, one 
subcell is occupied by the fiber material while the remaining three subcells are occupied by the 
matrix material. In most cases, the fiber material is at most transversely isotropic (as in the case 
of carbon fibers), while the matrix material is usually isotropic. In the form presented above, the 
original method of cells represents such a composite as a cubic material when in fact it is 
transversely isotropic. Aboudi (1987) overcame this problem by performing orientational 
averaging of the effective stiffness matrix in the x 2 - x 3 plane (see Figure 2.5) to yield a 
transversely isotropic effective stiffness matrix. This was done via the equation. 
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c-,.^)c‘{e)de, (2.113) 

0 

where c’(0) are the effective stiffness matrix components rotated by an angle 6 in the x 2 - x 3 
plane, and C* is the resulting transversely isotropic effective stiffness matrix. Brayshaw (1994) 

showed that the above averaging procedure results in an inconsistency in the stress concentration 
composition. That is, the geometrically-weighted subcell stresses do not sum to the cell global 
(cell) stresses as in equations (2.76) - (2.78) and (2.104) - (2.106). 

Brayshaw showed that the inconsistency is eliminated if orientational averaging is 
performed on stresses or strains rather than the effective stiffness terms. Further, it was 

determined that to ensure consistency, the integration limits should be — ^ to rather than 0 to 

n . In the reformulated original method of cells, Brayshaw averaging may be performed directly 
on the mixed concentration equations (2.70), (2.71), and (2.98) - (2.100). In the x 2 - x 3 plane, the 
rotation equations (which apply to both stresses and strains) are. 
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(2.114) 


(2.115) 


where c = cos 6 and s = sin 0 . Applying equation (2.1 14) to the subcell stresses yields, 
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The mixed concentration equations (2.70), (2.71), and (2.98) - (2.100) are then substituted into 
the above equations for the un-primed subcell stresses, yielding, for example, 
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The rotation equation (2.1 15) is then applied to the global strains, s # , thermal strains, Ar , 


and plastic strains, efj n , in the equations for each rotated subcell stress of the form above. This 
yields, for example, 
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Equations of the form above provide the rotated subcell stresses as a function of the global 
strains, thermal strains, and plastic strains, all of which are rotated as well. These rotated subcell 
stresses can now be averaged according to Brayshaw by use of the equation, 


n 
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This orientational averaging is performed for each subcell stress component to yield, 
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and the constants arising due to integration are. 
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The final step in the application of Brayshaw averaging involves simply replacing the 
appropriate terms in equations (2.70), (2.71), and (2.98) - (2.100) with the hatted terms above. 
Thus, the averaged global constitutive equations can be written as, 
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where the averaged effective stiffness components, coefficients of thermal expansion, and 
global plastic strain components are, 
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2.5 Classical Incremental Plasticity Theory 

The equations for calculation of the subcell plastic strain increments in the original 
method of cells subcells are here presented. The derivation follows that of Mendelson (1983) 
and simplifications provided by Williams and Pindera (1994) and Williams (1995). Omitting the 
designation {j3y) that identifies a given subcell in the original method of cells for notational 
simplicity, the subcell total strain components are given by the sum of the elastic strain 
components, the plastic strain components, and the plastic strain increments, 


e ij £ ij ■ 

The modified total strain components are defined as, 

* IJ * 

Combining equations (2.143) and (2.144) yields, 

£ ij ~ £ ij + . 

The mean dilatation is subtracted from (2.145) to give, 
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where, e‘ y and efj are deviatoric quantities defined as shown. Hooke’s elastic law is given by 


where, 
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is the stress deviator, and the Prandtl-Reuss equations are given by, 
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defj^a'ydX, (2.149) 

where, dX is the proportionality constant obtained from the consistency condition requiring that 
the stress vector remains on the yield surface during plastic loading. Eliminating the deviatoric 
stress using (2.147) and (2.149) yields, 
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and substituting for efj in equation (2.146) using (2.150) yields, 
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The equivalent modified total strain is defined as, 
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Substituting for in equation (2.152) using (2.151) yields. 
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where the definition of the effective plastic strain increment, de p ff has been 

Combining (2.151) and (2. 1 53) to eliminate the term l + — - — yields, 
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It can be shown that the proportionality constant, dk , can be expressed as. 
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is the effective stress. Substituting for dk in equation (2.153) using (2.156) results in, 
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where a modified proportionality constant, dX ' , has been defined. 

The plastic strain increments for the subcells in the original method of cells are calculated 
from equation (2.159). This equation represents a modification of the Prandtl-Reuss equations 
(2.149), such that the plastic strain increments are calculated from the modified total strain 
deviator rather than from the stress deviator. This form allows better convergence when it is 
employed in an iterative solution procedure. The terms in equation (2.159) are calculated from, 


e 



where, 



= H£ eff +&Y’ (2.160) 
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where E, H SP , and a Y are the elastic modulus, hardening slope, and yield stress for the material 
based on a bilinear stress-strain response. 


2.6 Solution Procedure in the Presence of Plasticity 

The objective of the analytical model is to predict the inelastic response of a material 
exhibiting periodic microstructure when subjected to global mechanical loading and a 
temperature change. This is done via the effective constitutive equation (2.40). In the present 
model, mechanical loading is imposed as global strain components, e y . The effective stiffness 

components, C*, and effective CTE components, or* , are determined from closed-form 
expressions involving the repeating unit cell geometry, the subcell compliance components, and 
the subcell CTEs. The compliance components and the CTEs of the heterogeneous subcells are 
determined using the original method of cells (see Section 2.3). One or more of the global strain 
components, ?. , are known from the imposed mechanical loading. The unimposed global strain 
components are calculated from the effective constitutive equation (2.40) with the use of global 
stress-free conditions. That is, for example, if e n is applied, typically a 0 =0 for ij * 1 1 . 
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When plasticity is present, the global constitutive equation (2.40) is nonlinear, and thus 
cannot be solved for the global stresses, a tJ , directly. The nonlinearity arises because the global 

plastic strain components, e? , depend on the subcell plastic strain components, s?j ' , which 

themselves depend implicitly on the subcell stress and strain components. The subcell stress 
and strain components can only be determined once the global constitutive equation is solved. 

Thus, to solve for the global stress components, cr y , iteration is necessary to find the correct z p 

for the imposed e . In addition, the mechanical loading in the form of imposed global strains, e , 
as well as any thermal loading, must be applied in an incremental manner. 

The three-dimensional subcell plastic strain components, , are calculated from the 

local subcell plastic strain components in the original method of cells, sf Pr \ At a given 
magnitude of the thermomechanical loading, these local plastic strains are expressed as the local 
plastic strains at the previous thermomechanical load plus an increment in the local plastic strains 
due to the increment in the thermomechanical load. This procedure was outlined by Mendelson 
(1983) and can be summarized as, 


~p(*) 

b ij 



previous 




m 


(2.162) 


where ds^ Pr ) are the increments in the subcell plastic strain components calculated using the 

original method of cells in conjunction with the classical incremental plasticity theory (Section 
2.5). Thus, once the local plastic strain increments are determined, the local plastic strains 
follow from equation (2.162), and the three-dimensional subcell plastic strains, upon which the 
unknown global plastic strains depend, can be calculated as well. 

The iterative procedure actually used in the model allows equation (2.40) to be bypassed 
during iteration. The mixed concentration equation (2.38) is used instead. The subcell and 
global stresses are calculated after convergence has occurred for a particular loading increment 
and are not active in the iteration procedure. The iterative procedure may be outlined as follows: 
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1) Apply a loading increment (i.e., a small increase in temperature or strain). 

2) Estimate the local plastic strains from equation (2. 1 62). 

3) Determine the three-dimensional subcell plastic strains using the original method of cells 
and determine the plastic terms in (2.38) from the equations in the appendix. 

4) Obtain the stresses in each three-dimensional subcell from equation (2.38). 

5) Determine three-dimensional subcell strains from the subcell constitutive equations (2.23). 

6) Apply these subcell strains to the original method of cells, from which new estimates of the 
local plastic strain increments are determined through the use of equation (2.159). 

7) Update the local plastic strains with equation (2.162) using the new local plastic strain 
increment values. 

8) Check for convergence of the local plastic strains. 

9) If convergence has been achieved, calculate the global stresses from equation (2.40) and go 
to step 1 . 

10) If convergence has not been achieved, go to step 3. 
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3. Modeling the Mechanical Response of 8H Satin C/Cu 

In this chapter, predictions of WCGMC (described in the previous chapter) are presented. 
In particular, the effects of unit cell microstructure and fiber volume fraction, porosity, residual 
stresses, and imperfect fiber-matrix bonding on the predicted monotonic tensile, compressive, 
and shear response of 8H satin C/Cu are examined. The effects are investigated in order to 
determine which are important and how each alters the model predictions. Those effects that are 
important will then be employed in Chapter 4 as the model predictions are compared with 
experiment. 

3.1 Effect of Unit Cell Microstructure and Fiber Volume Fraction 

One benefit of using WCGMC to model woven composites is the geometric flexibility 
offered by the embedded approach. As long as a repeating unit cell can be identified for a given 
heterogeneous material, the micro-scale geometry can be discretized into parallelepiped subcells 
and modeled with WCGMC. In addition, the averaged continuity conditions employed by the 
method of cells make precise geometric representation less important than it is for numerical 
finite-element or boundary-element models. It has been shown (Wilt, 1995) that for similar 
geometrical representations, the two-dimensional version of the generalized method of cells 
(GMC-2D) with 49 subcells matched elastoplastic finite-element predictions obtained using a 
1088 element mesh. For the above cases, the CPU time for the finite element model execution 
was 3550 times that required for GMC-2D execution. 

The geometry of the repeating unit cells for woven composites are often quite complex. 
The weave of fiber yams is inherently three-dimensional, and each different weave type has a 
different repeating unit cell (see Figure 1.3, for example). For an 8H satin weave, the repeating 
unit cell is shown (from above) in Figure 3.1. The unit cell is large, encompassing eight yams in 
each of the two directions, but it is the smallest rectangular repeating unit cell for this type of 
weave. A slightly smaller repeating unit cell can be identified if the rectangular shape 
requirement is relaxed, but for implementation in WCGMC, the unit cell must be a 
parallelepiped. 
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Figure 3.1. Schematic representation of the 8H satin weave and its repeating unit cell. 


For a composite reinforced with an 8H satin weave, the simplest geometric representation 
of the true three-dimensional unit cell is shown in Figure 3.2. This geometry is similar to the 
mosaic model of Chou and Ishikawa (1989) (see Figure 1.4) in the way it treats the fiber yam 
cross-over regions. The weave appears similar to a [0°/90°] laminate, with the 0° and 90° plies 
reversing stacking sequence at the yam cross-over points. Unlike the mosaic model, however, 
the geometry shown in Figure 3.2 represents the entire three-dimensional 8H satin unit cell, not 
just a two-dimensional section of it. This unit cell geometry will be referred to as the true mosaic 
model (TMM). Note that the TMM geometry incorporates pure layers of the copper matrix on 
the top and bottom of the unit cell, which mimics actual woven composites. The darkened top 
subcells in Figure 3.2 indicate the positions (in the plane of the weave) of the fiber yam cross- 
over regions. 

An additional point to consider when examining the unit cell geometry shown in Figure 
3.2 is the fact that WCGMC models the unit cell as a representative part of an infinite medium. 
That is, the unit cell repeats infinitely in each direction, not just in the plane of the weave. Thus, 
the geometry shown actually represents an infinite number of reinforcement weave layers, 
separated by regions of pure copper matrix. This contrasts with the models of Chou and 
Ishikawa (1989) and Naik and Ganesh (1992) which use lamination theory and thus model a 
plate, with free surfaces, reinforced with a single woven layer. 
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Figure 3.2. True Mosaic Model (TMM) unit cell geometry. 


The unit cell geometry shown in Figure 3.2 includes pure matrix layers above and below 
the weave, but no other regions of pure matrix. Thus, in order to vary the overall fiber volume 
fraction of the woven composite independently of the fiber volume fraction of the yam subcells, 
the thickness of the pure matrix layers is varied with respect to the weave subcell thickness. That 
is, the ratio x/h defined in Figure 3.2 is varied. Choice of the dimension w is arbitrary. In 
addition, the TMM geometry lacks fiber continuity at the yam cross-over points. In order to 
allow more freedom in varying the fiber volume fraction, as well as a more accurate 
representation of the yam cross-over geometry, a unit cell with considerably more subcells must 
be employed. Figure 3.3(a) and (b) show the next level of refinement for the 8H satin composite 
unit cell geometry. The unit cell geometry shown in Figure 3.3(a) is referred to as the continuous 
mosaic model (CMM), and the unit cell geometry shown in Figure 3.3(b), which explicitly 
includes rotated fibers in the cross-over regions, is called the directionally continuous mosaic 
model (DCMM). 
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Figure 3.3(a). Continuous Mosaic Model (CMM) unit cell geometry. 
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Figure 3.3(b). Directionally Continuous Mosaic Model (DCMM) unit cell geometry. 
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Clearly, the CMM and DCMM geometries provide superior representation of the woven 
composite compared to the TMM geometry. In addition, the CMM and DCMM geometries 
include pure matrix subcells within the reinforcement weave as well as those located in the top 
and bottom layers of the unit cell. This allows an additional geometric degree of freedom when 
tailoring the subcell dimensions to provide a particular overall fiber volume fraction. However, 
these refinements come at a computational cost. The TMM unit cell is comprised of 256 subcells 
(8x8x6), while the CMM and DCMM unit cells are each comprised of 1536 subcells (16x16x6). 
For implementation in WCGMC these numbers can be reduced to 192 (8x8x3) and 1280 
(16x16x5) by combining the top pure matrix layer of the unit cell with the bottom pure matrix 
layer of the unit cell above the unit cell being considered. WCGMC execution times are 
sensitive to the number of subcells comprising the unit cell considered. Indeed, prior to the 
reformulation of GMC-3D, the DCMM geometry could not be employed for elastoplastic 
simulations due to the large memory requirements and execution times required by the model 
(see Bednarcyk et al., 1997). However, the reformulation of the equations of GMC-3D has now 
made possible the use of the DCMM geometry for thermoinelastic modeling of woven metal 
matrix composites. 

For modeling 8H satin C/Cu, the fiber volume fraction of the fiber yam subcells was 
taken to be 65%. This value was determined via microscopic examination of actual C/Cu 
specimens and represents an average for many infiltrated yams in many composite specimens. 
Given the fiber volume fraction of the infiltrated fiber yams, the dimensions in the unit cells may 
be selected to obtain the desired overall fiber volume fraction. Table 3.1 provides the 
dimensions for the microstructures that will be considered. Note that results generated using the 
CMM geometry will not be presented since its microstructural accuracy is superseded by the 
DCMM geometry. In addition, macroscopic results generated using the CMM geometry are 
invariably nearly identical to those generated using the DCMM geometry (assuming the 
dimensions are the same). Two versions of the DCMM geometry will be considered, each with a 
different pure copper layer thickness. The DCMM #2 geometry is employed because of the 
upper bound of 43.3% on the overall fiber volume fraction using the DCMM #1 geometry. Both 
DCMM geometries model the actual microstructure of the 8H satin C/Cu composite specimens 
reasonably realistically. 
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V f (%) 

X 

■■ JWtKM 

a 

b 

TMM 

35 

0.0171 

0.02 

- 

- 


40 

0.0125 

0.02 

- 

- 


45 

0.0089 

0.02 

- 

- 


50 

0.0060 

0.02 

- 

- 

DCMM #1 

35 

0.01 

0.02 

i 

0.238 


40 

0.01 

0.02 

i 

0.083 

DCMM #2 

35 

0.005 

0.02 

i 

0.486 


40 

0.005 

0.02 

i 

0.300 


45 

0.005 

0.02 

i 

0.156 


50 

0.005 

0.02 

i 

0.040 


Table 3.1. Unit cell dimensions. 


The thermomechanical input properties required by the model for the carbon fiber and the 
copper matrix are given in Table 3.2. The properties of T-300 carbon fibers were used due to 
their similarity to VCX-1 1 fibers (Ellis, 1995). The fiber is treated as temperature-independent, 
elastic, and transversely isotropic. The copper matrix is treated as temperature-dependent, 
elastoplastic with linear hardening, and isotropic. 

3.1.1 Tensile Response 

Figures 3.4 - 3.7 show predicted monotonic tensile stress-strain curves for 8H satin C/Cu 
as represented by the different unit cell geometries for fiber volume fractions of 35%, 40%, 45%, 
and 50%. The composite effective elastic properties predicted by the model in each case are 
given in Table 3.3. Figures 3.4 - 3.7 indicate that the effect of unit cell geometry is minor for the 
predicted tensile response of the composite. The effect is most evident for the lowest fiber 
volume fraction, 35%. In this case, the TMM geometry provides the stiffest overall response, 
followed by the DCMM #1 geometry, and finally the DCMM #2 geometry. This trend is also 
apparent in the effective elastic moduli given in Table 3.3. In the elastic region of the composite 
stress-strain response, the TMM geometry is expected to provide the stiffest response. This is 
partly due to the fact that the TMM geometry has a larger number of yam subcells with fibers 
oriented in the loading direction than do the DCMM geometries. More importantly, at each fiber 
volume fraction, the TMM geometry has thicker pure copper layers than the more complex 
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T-300 Carbon Fibers 



e a 

(Msi) 

Et 

(Msi) 

Ga 

(Msi) 

va 


<*A 

(10-6/ o F) 

a T 

(10-6/op) 


29.42 

3.67 

6.40 

0.443 

0.05 

-0.389 

5.556 

Source 

1 

1 

1 

1 

1 

2 

2 


OFHC Copper 


Temp 

E 

V 

a 

°y 

H S p 

(°F) 

(Msi) 


(10"6/°F) 

(ksi) 

(Msi) 

70 

18.8 

0.35 

8.18 

10.3 

1.425 

100 

18.7 

0.35 

8.28 

9.97 

1.403 

300 

18.0 

0.35 

9.35 

9.37 

1.270 

500 

17.2 

0.36 

10.1 

8.95 

1.126 

700 

16.4 

0.37 

10.7 

8.42 

1.000 

900 

15.5 

0.375 

11.6 

6.15 

0.760 

1100 

14.5 

0.38 

12.0 

3.87 

0.521 

1300 

13.2 

0.38 

12.3 

2.25 

0.363 

1500 

11.7 

0.37 

12.6 

1.25 

0.282 

1700 

9.8 

0.35 

14.7 

0.27 

0.204 

Source 

3 

3 

3 

3/4 

3/4 


Table 3.2: Material properties for T-300 carbon fibers and OFHC copper. Sources: 1 - Derstine 
(1988); 2 - Naik and Deo (1992); 3 - Rocketdyne Materials Properties 
Manual (1987); 4 - NASA Lewis Research Center (1992). 


geometries. The stiffness of the pure copper is greater than the transverse stiffness of the carbon 
fibers (see Table 3.2), and since a large percentage of the fibers are oriented transverse to the 
loading direction, the thicker copper layers provide increased stiffness. The presence of the 
thicker copper layers also explains the stiffer response of the DCMM #1 geometry in the elastic 
region compared to the DCMM #2 geometry (see Table 3.2). 

Once yielding occurs in a given subcell, the “stiffness” of the copper in that subcell 
decreases dramatically. Thus, once yielding has occurred, the transverse stiffness of the elastic 
fibers is greater than the “stiffness” of the copper, and the pure copper layers become much less 
important. Now it is the pure copper subcells within the layers of the unit cell which represent 
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Figure 3.4. Predicted tensile stress-strain curves for 8H satin C/Cu, V f = 35%, 
using different unit cell geometries. 
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Figure 3.5. Predicted tensile stress-strain curves for 8H satin C/Cu, V f = 40%, 
using different unit cell geometries. 
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Figure 3.6. Predicted tensile stress-strain curves for 8H satin C/Cu, V f = 45%, 
using different unit cell geometries. 



Figure 3.7. Predicted tensile stress-strain curves for 8H satin C/Cu, V f = 50%, 
using different unit cell geometries. 
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Geometry 

v,(%) 

E (Msi) 

G (Msi) 

V 

TMM 

35 

16.04 

6.763 

0.2692 


40 

15.62 

6.734 

0.2544 


45 

15.19 

6.706 

0.2386 


50 

14.76 

6.677 

0.2215 

DCMM #1 

35 

15.69 

6.759 

0.2547 



15.47 

6.728 

* 0.2480 

DCMM #2 

35 

15.37 

6.759 

0.2413 


40 

15.12 

6.730 

0.2323 


45 

14.88 

6.700 

0.2242 


50 

14.66 

6.663 

0.2171 


Table 3.3 Effective elastic constants predicted by WCGMC for 8H satin C/Cu. 


the reinforcement weave that affect the behavior. After yielding, these pockets of pure copper 
serve as weak links, and bring about a more compliant post-yield response. The TMM geometry 
has no pure copper subcells within the weave and thus provides the stiffest post-yield response. 
Comparing the dimension b in Table 3.1 for the DCMM #1 and #2 geometries indicates that the 
pockets of pure copper are larger for the DCMM #2 geometry (see Figure 3.3b). Thus the 
DCMM #1 geometry has smaller weak links, and it provides a stiffer post-yield response. 

It is also apparent from Figures 3.4 - 3.7 that increasing the overall fiber volume fraction 
decreases the differences between the stress-strain curves predicted using the TMM and DCMM 
#2 geometries. As Figure 3.7 shows, these curves are essentially coincident for a fiber volume 
fraction of 50%. This too is explained by the size of the pure copper layers for the elastic region, 
and by the size of the weak link pure copper subcells for the post-yield region. The difference in 
the size of the pure copper layers between the TMM and DCMM #2 geometries decreases as the 
fiber volume fraction increases (see Table 3.1). Hence the tensile moduli predicted using the two 
geometries become closer in magnitude. Also, the size of the pure copper subcells within the 
reinforcement weave in the DCMM #2 geometry decreases with increasing fiber volume fraction 
(see Table 3.1). This explains the better agreement between the TMM and DCMM #2 
geometries at higher fiber volume fractions. 

Comparing Figures 3.4 and 3.5, however, indicates that, in the post-yield region, 
differences between the curves predicted using the DCMM #1 and DCMM #2 increase as the 
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fiber volume fraction is increased from 35% to 40%. This is because the difference between the 
b dimensions of the two geometries is much greater for a fiber volume fraction of 40% (see Table 
3.1). Hence, the difference in the size of the weak link copper subcells is greater for the 40% 
fiber volume fraction, which is apparent in the predicted stress-strain curves. Since the relative 
size of the pure copper layers does not change with fiber volume fraction for the DCMM #1 and 
DCMM #2 geometries (see Table 3.1), in the elastic region, increasing the fiber volume fraction 
has almost no effect on the difference between the predictions of the these geometries. 

3.1.2 Shear Response 

Figures 3.8 - 3.1 1 and Table 3.3 indicate that the unit cell geometry has a lesser effect on 
the predicted elastic shear response of the composite than on the predicted elastic tensile 
response. However, the effect of unit cell geometry is much more pronounced in the plastic 
region of the predicted shear stress-strain response than is the case in tension. The differences 
between the shear stress-strain curves predicted using the TMM and DCMM #2 geometries 



Figure 3.8. Predicted shear stress-strain curves for 8H satin C/Cu, V f = 35%, 
using different unit cell geometries. 
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Figure 3.9. Predicted shear stress-strain curves for 8H satin C/Cu, V f = 40%, 
using different unit cell geometries. 
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Figure 3.10. Predicted shear stress-strain curves for 8H satin C/Cu, V f = 45%, 
using different unit cell geometries. 
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Figure 3.11. Predicted shear stress-strain curves for 8H satin C/Cu, V f = 50%, 
using different unit cell geometries. 


decrease as the fiber volume fraction increases, as was the case in tension. Differences between 
the curves generated using the DCMM #1 and DCMM #2 geometries increase as the fiber 
volume fraction increases, which is also identical to the trend in the tensile predictions. Thus, the 
effects of the unit cell geometry on the predicted shear response of 8H satin C/Cu are the same as 
the effects on the predicted tensile response for the composite. The effects can also be explained 
exactly as they were for the tensile response, namely by the thickness of the pure copper layers 
for the elastic response, and by the size of the weak link pure copper subcells within the 
reinforcement weave for the post-yield behavior. 

In the elastic region, the shear modulus of pure copper is still greater than the axial shear 
modulus of the carbon fiber, but only slightly so (6.96 Msi vs. 6.40 Msi). Thus, in the elastic 
region, although the trend in stiffness based on geometry that was observed in tension is still 
present (see Table 3.3), the effect is much less apparent. Once yielding occurs, the shear 
modulus of copper drops as drastically as does the tensile modulus, however, the difference 
between the post-yield shear “stiffness” of copper and the shear modulus of the elastic fiber is 
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much greater than difference between the post-yield “stiffness” of the copper and the transverse 
tensile modulus of the fiber. Thus the magnitude of the effect of the unit cell geometry on the 
predicted composite post-yield behavior (caused by the weak link copper subcells) is much 
greater for shear loading than is the case for tension. As the figures and Table 3.3 indicate, 
however, the qualitative effects are identical. 

3.1.3 Summary of Fiber Volume Fraction Effects 

The effect of altering the overall fiber volume fraction on the predicted tensile and shear 
response of the composite using the DCMM #2 geometry is isolated in Figures 3.12 and 3.13. 
Note that in the elastic portion of the predicted stress-strain curves, increasing the fiber volume 
fraction leads to a more compliant response, while in the post-yield region, increasing the fiber 
volume fraction leads to a stiffer response. These trends are due to the low transverse tensile 
modulus and axial shear modulus of the carbon fibers (see Table 3.2). When the overall fiber 



Figure 3.11. Predicted tensile stress-strain curves for 8H satin C/Cu for different fiber volume 
fractions using the DCMM #2 unit cell geometry. 
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Figure 3.12. Predicted shear stress-strain curves for 8H satin C/Cu for different fiber volume 
fractions using the DCMM #2 unit cell geometry. 


volume fraction of the unit cell is increased, the volume fraction of the pure matrix subcells 
decreases, while the local fiber volume fraction of the infiltrated fiber yams remains constant. 
Since many fiber yams are oriented transversely to the loading direction in tension, and the 
transverse stiffness of the fibers is lower than that of pure copper, lowering the relative volume of 
pure copper lowers the tensile stiffness in the elastic region. This is true for shear loading as 
well, but since the axial shear modulus of the carbon fibers is only slightly lower than the shear 
modulus of copper, the effect is small. Once yielding of the copper matrix occurs, the transverse 
fibers are considerably stiffer in tension and shear than the plastically deforming copper. Thus, 
in the post-yield region of the stress strain-curves, lowering relative volume of pure copper (and 
thus increasing the overall fiber volume fraction) increases the stiffness. For reasons explained 
previously, the effect on the elastic response is greater in tension, while the effect on the post- 
yield response in greater in shear. 
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3.2 Effect of Porosity 

It was shown in the previous report on 8H satin C/Cu (Bednarcyk et al., 1997) that the 
composite plates fabricated for the investigation typically contained significant porosity. This 
porosity was not present in the pure copper matrix regions, but rather within the infiltrated fiber 
yams. Microscopic investigation of the composite yams revealed that this porosity typically 
occupied approximately 7% to 14% of the volume of the infiltrated fiber yams. 

In order to model its effect on the response of the woven composite, the porosity was 
treated via a reduction in the properties of the copper matrix within the infiltrated fiber yams. 
WCGMC was used to model the tensile response of pure copper with spherical void inclusions. 
This porous copper was modeled using a total of 343 subcells (7x7x7 cube). The overall volume 
fraction of the voids was set at 20% and then 40% so that a yam subcell with a fiber volume 
fraction of 65% infiltrated with these effective porous copper materials would have effective 
porosities of 7% and 14%, respectively. Figure 3.13 shows the predicted tensile response of the 
porous copper. The predicted response of copper with no voids and 1 0% porosity are plotted for 



Figure 3.13. Predicted tensile response for porous copper as a function of % porosity. 
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comparison. Clearly, porosity has a major effect on the predicted response of pure copper. Table 
3.4 provides a comparison of the effective room temperature properties for the porous copper 
determined by the model. 

To model the effect of porosity on the response of 8H satin C/Cu, the effective properties 
for the porous copper given in Table 3.4 can be used as input properties for the copper located 
within the infiltrated yam subcells in WCGMC. Figure 3.14 shows the predicted tensile 
response of 40% 8H satin C/Cu with 7% and 14% porosity within the infiltrated yam subcells. 


void content 

E(Msi) 

V 

G (Msi) 

Y (ksi) 

H SP (Msi) 

0% 

18.8 

0.350 

6.96 

10.3 

1.425 

10% 

14.3 

0.291 

5.55 

8.13 

1.11 

20% 

11.7 

0.260 

4.62 

6.61 

0.931 

40% 

7.26 

0.207 

3.01 

5.42 

0.597 


Table 3.4. Effective properties for porous copper determined by WCGMC. 
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Figure 3.14. Predicted tensile response for porous 8H satin C/Cu 
as a function of yam subcell % porosity. 
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The predicted response of the fully infiltrated composite is plotted for comparison. The DCMM 
#2 unit cell geometry was used for these predictions since it provides the most accurate 
geometric representation of the woven composite microstructure. From Figure 3.14 it is clear 
that porosity has a significant effect on the predicted tensile response of the woven composite. 
Higher levels of porosity lead to a more compliant overall predicted tensile response. The effect 
is noticeable even in the elastic region. Table 3.5 provides a numerical comparison of the 
effective tensile and shear properties of the porous composites as predicted by the model. While 
the predicted Poisson’s ratio for the composite changes only slightly as a function of the 
porosity, the predicted tensile modulus is reduced significantly by increasing fiber yam porosity. 

Figure 3.15 shows the predicted shear response for the porous woven composite. The 
predicted shear response for the fully infiltrated case is included for comparison. It is clear that 
porosity has an even greater effect on the shear response of the composite than on the tensile 
response. Once again, greater amounts of porosity lead to a more compliant predicted response 
for the composite. In addition, the predicted shear stress-strain curves in Figure 3.15 show that 
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Figure 3.15. Predicted shear response for porous 8H satin C/Cu 
as a function of yam subcell % porosity. 
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including yam porosity results in a much more gradual yielding of the composite compared to 
the fully infiltrated case. Table 3.5 provides a numerical comparison of the effective shear 
modulus for the porous composite. The data indicate that porosity, even when limited to the 
infiltrated fiber yams, significantly reduces the overall stiffness of 8H satin C/Cu. 


yam porosity 

E (Msi) 

G (Msi) 

V 

0% 

15.12 

6.730 

0.2323 

7% 

13.74 

6.155 

0.2170 

14% 

12.76 

5.541 

0.2110 


Table 3.5. Predicted effective properties for porous 40% 8H satin C/Cu. 


3.3 Effect of Residual Stresses 

In order to model the effect of residual stresses on the predicted response of 8H satin 
C/Cu, the cyclic capabilities of WCGMC were employed. Two loading cycles were 
implemented. The first is a cool-down from a temperature at which the composite is assumed to 
be globally and locally stress-free. This simulates the consolidation of the composite during 
production. A high stress-free temperature (SFT), close to the melting temperature of copper 
(1981 °F), would be chosen if the constitutive theory for copper modeled stress relaxation 
accurately. In order to use incremental plasticity, a lower effective SFT should be chosen to 
account for the relaxation of stresses in the copper matrix, which may occur at elevated 
temperatures during the cool-down. After the cool-down, the simulated mechanical loading is 
applied. During processing of the model predictions, the resulting mechanical stress-strain 
curves are shifted such that, at the starting point for the mechanical loading, the global strain on 
the composite is zero. This simulates the assumed zero strain starting condition for mechanical 
testing. Recall that the material properties of copper are temperature-dependent (see Table 3.2). 
WCGMC uses linear interpolation between input temperatures to determine the mechanical 
properties of the phases at the current temperature during thermal cyclic loading. 

As mentioned above, at elevated temperatures, copper is subject to viscoplastic behavior 
which can lead to relaxation of residual stresses. While the time-independent incremental 
plasticity constitutive theory is realistic for copper at room temperature, it has been shown that, 
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for the elevated temperatures associated with composite fabrication, a time-dependent plasticity 
theory should be employed (Bednarcyk and Pindera, 1996). One such viscoplasticity theory was 
developed by Freed et al. (1994) specifically for copper, commonly known as Freed- Walker 
viscoplasticity. While inclusion of Freed-Walker viscoplasticity in WCGMC is computationally 
prohibitive (due to the stiff nature of the equations which must be integrated), a version of the 
original method of cells embedded in classical lamination theory including this constitutive 
theory was developed by Bednarcyk and Pindera (1996). 

In order to determine a reasonable range of SFTs for 8H satin C/Cu for use in WCGMC, 
the aforementioned version of the inelastic lamination theory was employed to model the 
globally stress-free cool-down from 1700 °F of a cross-ply C/Cu laminate with a fiber volume 
fraction of 65% using the thermoelastic material properties given in Table 3.2. For the 
viscoplastic material parameters of copper which were employed see Bednarcyk and Pindera 
(1996). The cross-ply laminate provides the best simulation of an 8 H satin composite possible 
with lamination theory. A cool-down time of 900 seconds was employed to model the actual 
cool-down of 8H satin C/Cu composites during fabrication. Results, in the form of a plot of the 
effective stress in the matrix subcells in the plies vs. temperature are compared in Figure 3.16 
with predictions in which the incremental plasticity theory was used to model the constitutive 
behavior of the copper matrix. Note that, in this figure, the indices “12”, “21”, and “22” refer to 
the matrix subcell Py indices (see Figure 4.5). 

Clearly, utilizing plasticity in conjunction with an SFT of 1 700 °F vastly overpredicts the 
stresses in a unidirectional C/Cu composite. Results generated using Freed-Walker 
viscoplasticity indicate that a significant portion of the residual stresses predicted using 
incremental plasticity are not actually present at room temperature. Although the exact values of 
the residual stresses do not carry over from the cross-ply laminate to the woven composite, the 
qualitative conclusion does: an SFT of 1700 °F is unrealistic for WCGMC. 

In order to estimate a realistic value of the SFT for use in modeling 8H satin C/Cu, 
additional simulations of the cool-down of 65% cross-ply C/Cu laminate were performed. For 
these cases, the incremental plasticity constitutive theory was employed, and the SFT was varied. 
The predicted residual effective stress after the cool-down is plotted vs. the SFT in Figure 3.17. 
The SFT for which yielding does not occur during the cool-down is evident in this figure where 
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Figure 3.16. Effective stress in the subcells of a 65% volume fraction cross-ply C/Cu composite 
as modeled by the original method of cells embedded in lamination theory with the Freed-Walker 
viscoplasticity and incremental plasticity constitutive theories. 



Figure 3.17. Predicted residual stress in the subcells of a 65% cross-ply C/Cu laminate as a 
function of SFT, where incremental plasticity theory has been employed. 
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the slopes of the plotted curves change abruptly. This temperature is between 150 °F and 200 °F. 
Comparing Figure 3.17 with the Freed-Walker predictions in Figure 3.16, it is clear that in order 
to obtain realistic estimates of the residual stresses in C/Cu composites using incremental 
plasticity theory, an SFT lower than approximately 250 °F should be selected. This SFT is above 
the minimum SFT that causes yielding in the composite upon cool-down, and it will be used as 
an upper bound as the effect of varying the SFTs is investigated with WCGMC below. 

Predicted thermal expansion curves for 8H satin C/Cu resulting from cool-down from 
SFTs of 250 °F, 170 °F, and 120 °F are shown in Figure 3.18. Because of the low yield stress of 
the copper matrix and the high axial CTE mismatch between the fiber and the matrix (see Table 
3.2), the composite yields during simulated cooling from 250 °F. The composite remains elastic 
during simulated cooling from 170 °F and 120 °F. For an SFT of 250 °F, Figure 3.18 shows that 
by the time the composite reaches room temperature, the model predicts that it has undergone a 
significant global thermal contraction. During cooling, the copper matrix is placed in tension, 
while the carbon fibers are, for the most part, placed in compression. This is because the CTE of 
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Figure 3.18. Predicted thermal expansion curves for 40% 8H satin C/Cu for cool-downs from 
250 °F, 170 °F, and 120 °F, to room temperature. 



70 


the copper is much greater than the axial and transverse CTEs of the fibers (see Table 3.2). Thus, 
upon cooling, the copper matrix tends to contract more than the fibers, and the fibers constrain 
the contraction of the matrix. This results in local tensile stresses in the matrix and local 
compressive stresses in the fibers, while the total strains (i.e., thermal strains plus mechanical 
strains) in both phases are negative (i.e., contraction). 

The predicted effect of residual stresses due to cool-down from SFTs of 250 °F, 170 °F, 
and 120 °F on subsequent mechanical loading is illustrated in Figures 3.19 - 3.22. These figures 
compare the predicted tensile, compressive, and shear stress-strain curves for 8H satin C/Cu with 
a fiber volume fraction of 40% where different SFTs have been assumed. Note that an SFT of 70 
°F corresponds to the case with no residual stresses. This case is included in the figures for 
comparison. Once again, the DCMM #2 unit cell geometry has been employed for the same 
reason cited in the previous section. 
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Figure 3.19. Predicted tensile and compressive stress-strain curves for 40% 8H satin C/Cu with 
stress-free temperatures (SFTs) of 250 °F and 70 °F. 
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Figure 3.22. Predicted shear stress-strain curves for 40% 8H satin C/Cu with different stress-free 

temperatures (SFTs). 

Note that the predicted compressive response of the composite is explicitly investigated 
here for the first time. For cases without residual stresses, the model predictions for tension and 
compression are identical in magnitude, with opposite signs (i.e., the curves are inverted). When 
residual stresses are included, since the composite has non-zero internal stress and strain fields 
prior to mechanical loading, predictions for tension and compression are no longer identical. 
This is illustrated in Figure 3.19, which shows the predicted tensile and compressive stress-strain 
curves for SFT = 70 °F and SFT = 250 °F. Clearly, while the predictions for SFT = 70 °F are 
identical for tension and compression, the predictions for SFT = 250 °F are noticeably different. 
In particular, while yielding begins almost immediately in both tension and compression for SFT 
= 250 °F, the hardening is significantly greater in compression than in tension. 

Figures 3.20 - 3.22 indicate that increasing the SFT has two primary effects on the 
predicted response of 8H satin C/Cu: it leads to yielding at a lower stress and stiffer post-yield 
behavior. This is true for simulated tension, compression, and shear loading. Furthermore, the 
predictions are very sensitive to the choice of SFT, especially in tension. 
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In tension (see Figure 3.20), increasing the SFT from 70 °F to 120 °F decreases the 
predicted yield stress of the composite noticeably, and raises the post-yield slope slightly. The 
predicted yield stress is reduced further when the SFT is increased to 170°F. Recall that, upon 
cooling from 120 °F and 170 °F, predictions indicated that the composite remains elastic. Thus, 
upon initiation of simulated tensile loading, the composite still behaves elastically. For an SFT 
250 °F, however, the composite yielded during the simulated cool-down. Thus, upon initiation 
of simulated mechanical loading, regions of the composite are already “on” the yield surface, and 
the composite exhibits plastic deformation immediately. The result is that the predicted initial 
response of the composite is drastically changed for a relatively small change in the SFT. The 
effect of varying the SFT on the predicted yield stress for the lower SFT values is equally 
dramatic. It is also apparent from Figure 3.20 that increasing the SFT increases the predicted 
tensile post-yield slope for the composite, but only slightly so. This results in the apparent 
tendency of the predicted tensile stress-strain curves to converge at the higher strains. 

Qualitatively, the compressive predictions, shown in Figure 3.21, resemble the tensile 
predictions. Higher SFT causes lower predicted yield stresses and higher predicted post-yield 
slopes. However, the effect on the initial slope of the predicted stress-strain curve is 
considerably smaller in compression than in tension, while the effect on the post-yield behavior 
is greater in compression than in tension. Thus the predicted stress-strain curves cross in Figure 
3.21 while they do not cross in Figure 3.20. 

As in tension and compression, the model predictions for the shear stress-strain response 
of 8H satin C/Cu exhibit yielding at lower stresses and stiffer post-yield behavior as the SFT is 
increased (see Figure 3.22). Varying the SFT between 70 °F and 250 °F has almost no effect on 
the predicted initial shear response of the composite. However, as in tension, the predicted yield 
stress is significantly reduced as the SFT is increased. As in compression, due to the higher post- 
yield slopes associated with higher SFT, the predicted stress-strain curves cross, but they do so at 
much higher levels of stress and strain. 

The effects of varying the SFT discussed above can be explained by examining the 
evolution of plastic strain in the subcells of the 8H satin C/Cu repeating unit cell during 
simulated mechanical loading. Figure 3.23 is a plot of the predicted average effective plastic 
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Figure 3.23. Predicted average effective plastic strain in several three-dimensional subcells as a 
function of the applied global strain for SFTs of 250 °F and 70°F. The numbers on the plot 
correspond to the subcells numbered in the diagram above. 


strain in several subcells vs. the applied global strain for simulated tensile loading for SFTs of 70 
°F and 250 °F. The diagram at the top of the figure identifies the subcells to which the curves in 
the figure correspond. The particular subcells chosen, though all from the same cross-section of 
the unit cell, are representative of most subcells in the repeating unit cell. 

The reasons for the earlier yielding and stiffer post-yield behavior predicted in tension as 
the SFT was increased are readily apparent in Figure 3.23. For an SFT of 250 °F, many subcells 
have yielded during the cool-down and begin to accumulate additional plastic strain immediately 
upon application of the simulated tensile loading. This contrasts with the curves that represent 
an SFT of 70 °F which indicate that the subcells remains elastic (with zero plastic strain) 
initially. As the SFT is increased from 70 °F, even if no subcells yield upon cool-down, residual 
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stresses in the subcells become greater. Thus, upon application of the simulated tensile loading, 
certain subcells yield sooner, lowering the composite yield stress. This explanation is valid for 
predicted compressive and shear results, which exhibited the identical trend in yield stress based 
on SFT, as well. 

The greater post-yield slopes which resulted from increased SFTs can also be explained 
via Figure 3.23. Certain subcells in the repeating unit cell (#1, #3, and #1 1) remain elastic to a 
higher level of the applied global strain for an SFT of 250 °F compared to an SFT of 70 °F. This 
is because, after the simulated cool-down, the state of residual stress in the matrix of a given 
subcell may be unfavorable for accumulation of plastic strain in response to the global loading. 
That is, the global tensile loading may cause the state of stress in the subcell to pass through a 
portion of the yield surface, and load elastically. This is, in effect, due to the highly non-uniform 
state of residual stress and strain in the repeating unit cell that results from the simulated cool- 
down, allowing a wide range of stress states in the subcells. Since some subcells load elastically 
longer and accumulate plastic strain at a lower rate, higher SFT gives rise to stiffer post-yield 
behavior. The reason for the increase in the post-yield stiffness as the SFT is increased for 
compression and shear is identical. 

Comparing the trends based on choice of SFT for tension and compression, it was 
observed that the effect of SFT on yield stress is greater for tension than compression, while the 
effect of SFT on the post-yield behavior was greater in compression than tension. Recall that, 
after the simulated cool-down, the matrix in most subcells is in tension. Thus, upon application 
of simulated global compression, many more subcells load elastically initially compared to 
simulated global tension. Hence, for higher SFTs, the initial slope of the predicted compressive 
stress-strain curve is not decreased significantly from the initial elastic slope of the predicted 
compressive stress-strain curve for SFT = 70 °F. The curves in Figure 3.21 therefore do not 
initially separate as much as those in Figure 3.20. However, since so many subcells do initially 
load elastically for high SFT in compression, the stiffening effect on the post-yield behavior is 
increased compared to tension. 

The modeled effect of SFT on the shear response falls between the effects of SFT on 
tension and compression. In shear, the effect of SFT on the predicted yield stress is greater than 
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the effect in compression, but not as dramatic as the effect in tension. Conversely, the effect of 
SFT on the post-yield behavior is greater than the effect in tension, but smaller than the effect in 
compression. This is because, since no in-plane shear stress arises in the composite during 
simulated cool-down, the state of residual stress is not biased against plastic strain accumulation, 
as was the case in compression. Thus, the large number of subcells which initially load 
elastically that affect the predicted compressive response of the composite are absent during 
simulated shear loading. However, upon application of the simulated shear loading, plastic strain 
initially accumulates more slowly than in tension, also because the state of residual stress 
involves no in-plane shear. This slow plastic strain accumulation rate causes the predicted shear 
stress-strain curves in Figure 3.22 to nearly coincide initially, even for the SFT of 250 °F. 

3.4 Effect of Imperfect Fiber-Matrix Bonding 

Imperfect fiber-matrix bonding is incorporated in WCGMC on the local level, through 
the original method of cells. Hence, the imperfect bonding that is modeled occurs between the 
fibers and the matrix within the infiltrated fiber yams. As described in Section 2.3, the imperfect 
fiber-matrix bonding is accounted for by a normal debonding parameter, R n , and a tangential 
debonding parameter, R , . These parameters represent an effective compliance, in the two 
directions, of the interface between the fiber and matrix in the original method of cells. 
Debonding parameter values of zero correspond to a perfectly bonded interface, while as the 
debonding parameters become large, total interfacial debonding is approached. The parameters 
may be varied independently as well. For instance, if R n = 0, and R, is large, the interface is 
perfectly bonded in the normal direction, but completely debonded tangentially. This condition 
corresponds to perfect sliding at the fiber-matrix interface (Aboudi, 1988). 

3.4.1 Elastic Moduli 

The effect of imperfect fiber-matrix bonding on the predicted elastic properties for 8H 
satin C/Cu with a fiber volume fraction of 40% will be examined first. As before, the DCMM #2 
unit cell geometry will be utilized for this purpose. Figures 3.24, 3.25, and 3.26 show the effects 
of varying the debonding parameters on the predicted tensile modulus, shear modulus, and 
Poisson’s ratio of 8H satin C/Cu. In Figures 3.24 and 3.26, each curve corresponds to a constant 
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Figure 3.26. Predicted effective Poisson’s ratio of 40% 8H satin C/Cu for different values of the 

debonding parameters R n and^, . 


value of R, , given in the legend, while the value of R n is varied along each curve. In figure 
3.25, each curve corresponds to a constant value of R n , given in the legend, while the value of 
R, is varied along each curve. Figure 3.24 indicates that increasing either debonding parameter 
results in a significantly lower effective tensile modulus of the composite. The most drastic 
reduction in tensile modulus occurs when the debonding parameters are small, while as each 
debonding parameter becomes large, the tensile modulus converges to a constant value. For the 
fully debonded cases (i.e., R n — > oo, R t — > oo) (not plotted in Figure 3.24), the predicted tensile 
modulus for 40% 8H satin C/Cu is 1 1.12 Msi, which represents a reduction of 26.4% from the 
perfectly bonded case. Clearly, imperfect fiber-matrix debonding can have a significant impact 
on the predicted tensile modulus for 8H satin C/Cu. Figure 3.25 indicates that, while increasing 
R, results in a significantly lower predicted shear modulus for 8H satin C/Cu, the normal 
debonding parameter, R n , has no effect. 
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Figure 3.26 shows that the effect of increasing the normal debonding parameter on the 
predicted Poisson’s ratio for 8H satin C/Cu is similar to the effect on the predicted tensile 
modulus. Namely, the Poisson’s ratio is greatly reduced initially by increasing R n and then 
levels off as R„ continues to increase. The effect of increasing R, on the predicted Poisson’s 
ratio is initially similar to the effect of increasing R, on the tensile modulus as well. However, 
the predicted Poisson’s ratio appears to converge to a constant value at a lower value of R, than 
does the predicted tensile modulus. In addition, the predicted Poisson’s ratio for R„ =5x 10' 6 
first decreases, but then increases slightly as R, is increased from 5xl0' 7 to 5 x 1 CT 6 . This 
appears to be due to the fact that convergence of the predicted Poisson’s ratio occurs at a lower 
value of R n as well, thus the curves in Figure 3.26 are, in effect, shifted to the left for the higher 
R , values, compared to the curves in Figure 3.24. 

The effect of R, on the normal behavior of the composite (including the effective tensile 
modulus, shown in Figure 3.24) is mainly due to the averaging procedure employed in the 
original method of cells (see Section 2.4). If Brayshaw averaging were not performed in the 
original method of cells, the effects of the normal and tangential debonding parameters would be 
completely decoupled for unidirectional composites. That is, the normal debonding parameter 
would affect only the normal response of the unidirectional composite, while the tangential 
debonding parameter would affect only the shear response. In the presence of rotated fibers, 
however, the tangential and normal debonding effects will be linked when averaging is not 
performed, but only in the plane of fiber rotation. That is, for fibers rotated in the x, - x j plane, 
the normal and tangential debonding effects would be linked between the normal stresses a„ and 
<?jj (no summations), and the shear stress <j tJ . However, in the case of 8H satin C/Cu as 
represented by the DCMM #2 geometry (see Figure 3.2(b)), all rotated fibers present are rotated 
in the x ] -x 2 and -* 3 planes. No fiber rotation is present in the plane of the weave, the x 2 -x 3 
plane. Thus, for 8H satin C/Cu, without Brayshaw averaging, the effects of the normal and 
tangential debonding would be expected to be decoupled for simulated in-plane shear loading. 
For in-plane normal loading, since small out of plane shear stresses arise in some subcells, the 
tangential debonding parameter would be expected to have a small effect. 
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Brayshaw averaging, in effect, links the debonding parameter effects for the in-plane 
normal response and in-plane shear response. This can be seen clearly in equations (2.130) - 
(2.133) and (2.137). The averaged terms which apply to the in-plane normal behavior in the 

original method of cells, B ( / y) , X ( / x) , X^”\ B Bf and Xf 1 ^ , 

are related to the un-averaged in-plane shear terms, A and , which are affected by the 
tangential debonding parameter (see equations (2.104) and (2.98)). Hence, the normal response 
of the composite (and thus the effective tensile modulus) is affected by R , . The effect of R, on 
the normal response of 8H satin C/Cu is thus an artifact of the averaging procedure, rather than a 
true effect of the sliding interface modeled by R , . However, since it affects the normal response 
in the same way as does R n , the utility of the compliant interface imperfect bonding model 
should not be diminished by this fact. 

It was observed in Figure 3.25 that R n does not affect the predicted in-plane shear 
modulus of the composite, even though Brayshaw averaging is employed. The reason for this 
can be determined by once again examining the averaging equations. Equations (2.137), (2.73), 
(2.74), (2.68), (2.69), and (2.70) show that, when Brayshaw averaging is performed, the averaged 

in-plane shear terms, and A /W ^ r ’^ , are influenced by the un-averaged normal terms, which 
are affected by the normal debonding parameter, R n . However, the effect of R n is identical for 
the in-plane normal terms in the 2 and 3 directions. Examining equations (2.137) indicates that, 
when this is the case, the effect of the normal debonding parameter on the averaged in-plane 

shear terms, A^ and A p '^ r ^ 9 will subtract out due to the form of the equation. Thus, as was 
observed in Figure 3.25, R n does not affect the predicted in-plane shear response, which includes 
the effective in-plane shear modulus, of the composite. 

3.4.2 Mechanical Response 

The effect of imperfect fiber-matrix bonding on the predicted tensile response of 40% 8H 
satin C/Cu is shown in Figures 3.27, 3.28, and 3.29. In these figures, each curve represents the 
predicted stress-strain response for a different value of R n , given in the legend. In each of the 
three figures, a different value of R, was used; 0, lxlO' 7 , or 1. Using a value of 1 for either 
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Figure 3.27. Predicted tensile stress-strain curves for 40% 8H satin C/Cu for different values of 

R n with R, =0. 
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Figure 3.28. Predicted tensile stress-strain curves for 40% 8H satin C/Cu for different values of 

R„ with R= lxlO" 7 . 
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Figure 3.29. Predicted tensile stress-strain curves for 40% 8H satin C/Cu for different values of 

R n with R, = 1 . 


debonding parameter is large enough to correspond to the completely debonded case. These 
figures indicate that the effect of increasing either debonding parameter on the predicted tensile 
response of the composite is to significantly lower the elastic stiffness, yield point, and post-yield 
stiffness. 

Figure 3.30 shows the effect of imperfect fiber-matrix bonding on the predicted shear 
response of 8H satin C/Cu. An R„ value of zero was used in this case, and R, was varied to 
produce the three curves. As discussed earlier, R n does not affect the predicted in-plane shear 
response of the composite. As the figure indicates, the effect of varying R, on the predicted 
shear response of the composite is substantial; even greater than is the effect of varying R n on 
the tensile response of the composite. In addition to the significant reduction of the predicted 
shear modulus, yield point, and post-yield stiffness, increasing the tangential debonding 
parameter smoothes out the yielding of the composite, making it less abrupt. The effect of R, on 
the predicted shear response of the composite is greater than the effect of R n on the predicted 
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Figure 3.30. Predicted shear stress-strain curves for 40% 8H satin C/Cu for different values of 

R , with R n = 0. 

tensile response of the composite because, while R n only directly affects one half of the fiber 
yams (those oriented transverse to the loading directions), R, affects all of the fiber yams. 

3.5 Summary of Modeled Effects 

Unit Cell Microstructure: Utilizing the more refined DCMM unit cell geometry gave rise to a 
slightly more compliant overall predicted tensile response for 8H satin C/Cu compared to results 
generated using the TMM geometry. The effect of unit cell refinement on the predicted elastic 
shear response of the composite was similarly small, but the effect was much more significant on 
the predicted shear post-yield response. This is due to the large difference between the axial 
shear modulus of the carbon fibers and the post-yield shear “stiffness” of the copper matrix. 


Fiber Volume Fraction: Increasing the overall fiber volume fraction of the composite gave rise 
to a more compliant predicted elastic response, but a stiffer post-yield response in both tension 
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and shear. This is due to the drop in the effective “stiffness” of the copper associated with 
yielding. Once again, the effect on the post-yield behavior was greater for shear than tension. 

Porosity: Increasing the porosity within the infiltrated fiber yams of 8H satin C/C u gave rise to a 
significantly more compliant overall predicted tensile response. The effect was qualitatively 
similar, but of larger magnitude, for the predicted shear response of the composite. 

Residual Stresses: Through utilization of a previously developed version of the original method 
of cells embedded in lamination theory utilizing a viscoplastic theory, it was shown that for high 
stress-free temperature (SFT) values, the model (based on incremental plasticity constitutive 
theory) significantly overpredicts the residual stresses in the composite. In order to obtain 
realistic residual stresses prior to simulated mechanical loading, it is necessary to use an SFT 
below approximately 250 °F. Small changes in the SFT chosen for the composite had a 
significant effect on the predicted mechanical response of 8H satin C/Cu. Higher SFT values 
lead to lower predicted yield stresses for the composite, but stiffer post-yield behavior. The 
effect on yield stress was greatest for simulated tensile loading and smallest for simulated 
compressive loading, while the effect on the post-yield behavior was greatest for simulated 
compressive loading and smallest for simulated tensile loading. For simulated shear loading, the 
magnitudes of both effects were intermediate. 

Imperfect Fiber-Matrix Bonding: Increasing either the normal or the tangential debonding 

parameter produced a significantly more compliant overall predicted tensile response for 8H 
satin C/Cu. Increasing the tangential debonding parameter produced a significantly more 
compliant overall predicted shear response for the composite, while the normal debonding 
parameter had no effect on the predicted shear response. The effect of the tangential debonding 
parameter on the predicted shear response was greater than the effect of either debonding 
parameter on the tensile response. 



4. Model-Experiment Correlation 


In this chapter, the mechanical test results for 8H satin C/Cu reported by Bednarcyk et al. 
(1997) are summarized, and the predictions of the developed micromechanics model are 
compared with the experimental data. Comparison is performed for monotonic and cyclic 
tension, compression, and shear tests, as well as combined tension-compression and combined 
compression-tension tests. By comparing the model to experimental data, the effectiveness and 
accuracy of the model can be evaluated. Furthermore, by seeking the best possible agreement 
with the model, additional insight can be gained on how and to what extent the previously 
identified factors affect the actual behavior of the composite. 

4.1 Summary of Experimental Results 

For a complete presentation and discussion of the mechanical tests and results for 8H 
satin C/Cu see Bednarcyk et al. (1997). Typical monotonic tensile and compressive stress-strain 
curves for the composite are plotted in Figure 4.1 along with the stress-strain response of the 
carbon fiber and the copper matrix. Note that, in tension, the composite is more compliant than 
pure copper and yields at a lower stress as well. The composite response exhibits stiffening at 
moderate strains and failed via breakage of the longitudinal fibers. In compression, the stress- 
strain response of the composite is quite different. The initial compressive stiffness is somewhat 
higher than the tensile stiffness, and the onset of yielding is more gradual in compression. 
Compressive failure occurred via micro-buckling of the layers of the reinforcement weave, 
resulting in a lower strength and strain to failure compared to tension. It should be noted that a 
specialized compressive test fixture was used for the compressive (and combined) tests. Hence, 
artifacts caused by this compressive fixture may be present. 

Since the fiber-matrix bonding in C/Cu composites is poor, small amounts of titanium or 
chromium were added to the matrix of some 8H satin C/Cu plates fabricated for the experimental 
investigation. These alloying additions improve the wetting between the copper and the carbon 
fiber and thus improve the bonding (Ellis, 1992; DeVincent, 1995). Previous work on C/Cu 
composites indicated that the C/Cu - 0.7 wt. % Cr (C/Cu-Cr) composite would have the best 
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Figure 4. 1 . Typical monotonic tension and monotonic compression response of 

8H satin C/Cu-alloy composites. 


fiber-matrix bonding, followed by the C/Cu - 0.5 wt. % Ti (C/Cu-Ti) composite, and finally by 
the pure copper matrix (C/Cu) composite. By examining differences in the responses of the 8H 
satin C/Cu composites with the three matrix alloy types, the effects of poor fiber-matrix bonding 
can thus be examined. 

The effect of matrix alloy type (and thus the effect of poor fiber-matrix bonding) on the 
tensile and compressive response of 8H satin C/Cu is shown in Figure 4.2. In compression, 
improved fiber-matrix bonding gives rise to a stiffer overall response. In tension, however, this 
trend is reversed. The tensile response of the C/Cu-Cr composite, with its superior fiber-matrix 
bonding, is most compliant. 

The in-plane shear response of 8H satin C/Cu with the three different matrix alloy types 
is shown in Figure 4.3. The shear test were of the Iosipescu type, employing a doubly-notched 
specimen and a specialized test fixture (see Bednarcyk et al. (1997) for details). Shear yielding 
occurred at low stresses, and a great deal of inelastic strain accumulated during the tests. Shear 
tests were performed in displacement control and were stopped when contact of components of 
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1.2. Typical monotonic tension and monotonic compression stress-strain curves for 8H 

satin C/Cu-alloy specimens. 
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Figure 4.3. Typical monotonic shear stress-strain curves for 8H satin 

C/Cu-alloy specimens. 
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the Iosipescu test fixture was imminent. Hence, failure in shear was not achieved. The trend 
based on matrix alloy type in shear is identical to that in compression, with improved fiber- 
matrix bonding giving rise to a stiffer overall response. 

Typical cyclic tensile, compressive, and shear test results are shown in Figures 4.4, 4.5, 
and 4.6. The hysteresis loops that characterize the cyclic tensile and compressive response of the 
composite at higher stress levels are caused by a combination of frictional sliding at the fiber- 
matrix interface and kinematic hardening rather than damage. The hysteresis loops present in the 
cyclic shear stress-strain curve (Figure 4.6) are considerably smaller. 

Figures 4.7 and 4.8 show the results of typical combined tension-compression (T-C) and 
combined compression-tension (C-T) tests on 8H satin C/Cu. Note that, in both combined tests, 
a great deal of permanent tensile strain was present after unloading in tension. Upon unloading 
from compression, no permanent compressive strain was evident. 
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Figure 4.4. Typical cyclic tensile stress-strain curve for 8H satin C/Cu. 




90 


2° j 



Figure 4.7. Typical combined tension-compression test results for 8H satin C/Cu. 
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Figure 4.8. Typical combined compression-tension test results for 8H satin C/Cu. 


91 


4.2 Monotonic Correlation 

In order to realistically model the tensile, compressive, and shear stress-strain response of 
8H satin C/Cu, it was desirable to use a well-refined geometric representation of the composite 
unit cell. Thus the DCMM #2 geometry was employed. An overall fiber volume fraction of 
40% was used for the simulations since this value is an approximate average for all 8H satin 
C/Cu specimens, while an infiltrated yam fiber volume fraction of 65% was used. As was shown 
in Section 3.1, small variations in fiber volume fraction should not have a significant effect on 
the predicted elastic response of the composite. Furthermore, the effect of slight fiber volume 
fraction variations (on the order of 1 % or 2 %) on the predicted post-yield behavior of the 
composite are insignificant, compared to the scatter in the experimental results, as well. 
Likewise, a yam porosity of 14% was employed throughout the correlation simulations since this 
value is realistic and representative. 

For simulations of tests involving the compressive gripping apparatus (monotonic and 
cyclic compressive tests and combined tests), a constraint on the lateral deformation of the 
composite specimen was imposed. The test fixture used for these tests requires a specialized 
non-dogboned specimen. The specimen is a rectangular plate measuring 1 inch by 4 inches, but 
1 inch of either end is held rigidly in the gripping apparatus. This leaves only a 1 inch by 2 
inches test section outside of the grips. It is suspected that, since the specimen is not dogboned 
and the aspect ratio of the test section is small, the grips constrain the 8H satin C/Cu specimen 
from expanding naturally in the in-plane direction transverse to the compressive loading. 
WCGMC was used to model this grip constraint effect, when appropriate, by including a global 
loading condition which, upon in-plane tensile or compressive normal loading, requires the 
global mechanical strain component transverse to the loading direction to be zero, rather than the 
stress component in the transverse direction (as is normally done). The effect of including this 
grip constraint effect on the predicted monotonic compressive response of the composite is 
shown in Figure 4.9. The constraint causes a slightly stiffer overall predicted response. 

With the above model input parameters, as well as the temperature dependent fiber and 
matrix material properties (see Table 3.2) fixed, three parameters remained to be chosen: the 
debonding parameters, R n and R , , and the stress-free temperature. As shown in Sections 3.3 
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Figure 4.9. Predicted compressive stress-strain curves for 40% 8H satin C/Cu with and without 

the constraint effect from specimen gripping. 


and 3.4, these parameters have a major impact on the model predictions, and thus must be chosen 
carefully. In choosing realistic values of R n and R , , the nature of the imperfect bonding model 
must be considered. R„ and R, represent the normal and tangential compliances of an effective 
interfacial layer (see Section 2.3) in the infiltrated fiber yams of the woven composite. Thus, 
non-zero values of these parameters give rise to jumps in the displacement components at the 
fiber-matrix interface which are proportional to the normal and shear stress components at the 
interface. This compliant interface model is quite realistic for interfacial shear stresses. If a 
positive or negative shear stress arises at an imperfectly bonded interface, a proportional jump in 
a tangential displacement component results (governed by /?,). Thus, utilization of R, for 
simulation of cyclic thermo-mechanical tests on 8H satin C/Cu is warranted. 

Unfortunately, the compliant interface model is somewhat unrealistic for interfacial 
normal stresses. If the interfacial normal stress is tensile, the model is representative, giving rise 
to a positive jump in the displacement component normal to the interface. However, if the 
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interfacial normal stress is compressive, the compliant interface model breaks down. It now 
predicts a negative jump in the normal displacement component at the poorly bonded interface, 
when in fact, the interface should simply remain closed, with continuity of the normal 
displacement component. Clearly, utilization of R n is unrealistic for simulated loading in which 
fiber-matrix interfaces are placed in compression. In the case of 8H satin C/Cu as represented by 
the DCMM #2 geometry, almost any simulated loading places some fiber-matrix interfaces in 
compression. Furthermore, upon application of a simulated cool-down, the fiber-matrix 
interfaces are placed in a state of residual compression since the CTE of copper is greater than 
the transverse CTE of the fiber. Thus, even when the interface is placed in tension during 
subsequent simulated mechanical loading, the interface should remain closed initially. In order 
to keep the simulations as realistic as possible, an R n value of zero was used for all simulations, 
disallowing the normal discontinuity in the displacement at the interfaces. 

Recall from Section 3.4 that, due to the averaging procedure performed in the original 
method of cells, R, does affect the normal response of the infiltrated fiber yams. Thus, inclusion 
of R, will have a noticeable effect on the tensile and compressive simulations. Furthermore, the 
need for imperfect bonding in order to obtain reasonable model-experiment correlation in shear 
was demonstrated by Bednarcyk et al. (1997). In this report, preliminary shear predictions (made 
by WCGMC) for perfectly bonding 8H satin C/Cu were shown to be in poor agreement with the 
experimental results. In all cases, the model significantly overpredicted (by almost a factor of 
two) the shear stress-strain response of the composite, while model predictions for tension and 
compression were reasonable. Hence, the motivation for employing R t , but neglecting R n , is 
clear. Since interfacial bonding is known to vary between the matrix alloy types, a different R , 
value was employed for each. R, values of 2xl0" 7 , 4x1 0‘ 7 , and 6x1 O' 7 were used for the C/Cu- 
Cr, C/Cu-Ti, and C/Cu composites, respectively. These values were chosen to provide 
reasonable agreement with experimental shear results rather than experimental tensile or 
compressive results because R , affects the shear response of the composite directly and most 
significantly. 

The remaining parameter to be chosen is the composite’s stress-free temperature. As 
shown in Section 3.3, utilization of a high stress-free temperature is unrealistic. Utilization of a 
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single stress-free temperature for all three matrix alloy types is unrealistic as well. It is likely 
that, due to the different levels of imperfect fiber-matrix bonding, each matrix alloy type results 
in a different stress-free temperature for the composite. Since the fiber-matrix bonding in the 
C/Cu composite is inferior to that in the C/Cu-Cr composite, for example, residual stresses would 
tend to relax to a larger extent since less restraint would be imposed by the fiber on the matrix. 
This greater relaxation cannot be modeled well by the compliant interface debonding model 
because it is caused predominantly by longitudinal displacement discontinuities between the 
fiber and the matrix. That is, a jump occurs in the longitudinal (fiber direction) displacement 
component at the poorly-bonded interface within the infiltrated fiber yams. This is because the 
large CTE mismatch between the fiber and the matrix in C/Cu occurs in the longitudinal 
direction (8. 18x1 O'VF vs. -0.389x 1 O' 6 /®?) rather than the transverse direction (8.18xlO' 6 /°F vs. 
5.556x1 0' 6 /°F). The original method of cells treats the longitudinal displacement of the entire 
infiltrated fiber yam as constant, and the imperfect bonding model offers no jump in the 
longitudinal displacement. Indeed, it was observed that choice of the debonding parameters had 
little effect on the simulated cool-down for a given stress-free temperature. 

Thus, it appears that the imperfect interfacial bonding model employed is insufficient to 
realistically model the effect of imperfect bonding on residual stresses, necessitating the use of a 
different stress-free temperature for each matrix alloy type. Values of stress-free temperature 
were chosen to be realistic (as discussed in Section 3.3), to provide good correlation with 
experiment, and to reflect the proper trend based on the degree of interfacial bonding for each 
matrix alloy type (i.e., C/Cu-Cr has the highest stress-free temperature, while C/Cu has the 
lowest). Stress-free temperatures of 170 °F, 120 °F, and 70 °F were utilized for C/Cu-Cr, C/Cu- 
Ti, C/Cu, respectively. It is thus assumed here that the C/Cu composite is completely stress-free 
at room temperature prior to mechanical loading. Though residual stresses in the C/Cu 
composite are probably quite low, assuming that they are zero is a simplification. 

Figure 4.10 compares the predicted tensile stress-strain response of 8H satin C/Cu-alloy 
composites, where the parameters presented above have been used, with three experimental 
tensile stress-strain curves for composites with each matrix alloy type. Clearly, the overall 
correlation between model and experiment is acceptable. In the elastic region, the three 
predicted stress-strain curves are nearly identical, and they agree well with the experimental 
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Figure 4.10. Comparison of model predictions (V f = 40%) and experimental data for the 
monotonic tensile stress-strain response of 8H satin C/Cu-alloy composites. 


curves. At the onset of yielding, the trend based on matrix alloy type is correctly predicted. The 
simulated 8H satin C/Cu-Cr composite ( R , = 2 x 10~ 7 , SFT =170 °F), yields first, followed by 
the simulated 8H satin C/Cu-Ti composite (R, = 4xl0" 7 , SFT =120 °F), and finally by the 
simulated 8H satin C/Cu composite ( R, = 6xl0" 7 , SFT = 70 °F). This correctly predicted trend is 
brought about by the stress-free temperatures chosen for the composites. As shown in Section 
3.4, the trend in R, chosen for the composites would cause the opposite trend in yielding. 
Hence, the effect of the poor interfacial bonding on the residual stresses in the composite (as 
modeled by choice of stress-free temperature) dominates the trend based on matrix alloy type at 
low strain. The direct effect of poor interfacial bonding on the predicted tensile response of the 
composite is clearly overshadowed. 

Examining the hardening behavior exhibited by the predicted tensile stress-strain curves 
in Figure 4.10, it clear that the model predicts the opposite trend based on matrix alloy type 
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compared to the experimental results. As was shown in Section 3.3, employing a high SFT in 
the model gives rise to a highly non-uniform residual stress states in the composite. Thus, during 
the subsequent simulated tensile loading, a number of subcells have stress states that are initially 
unfavorable for accumulation of additional plastic strain. The result is stiffer predicted post-yield 
behavior as the SFT is increased. This effect works in concert with the direct effect of decreased 
R t values, and as shown in Section 3.4, causes stiffer predicted post-yield behavior as well. 
Thus, while the experimental tensile stress-strain curves for each matrix alloy type continue to 
diverge at higher strains, the predicted curves converge. In order to compensate for this 
discrepancy, a more realistic imperfect bonding model would be required. 

Figure 4.1 1 compares the predicted compressive stress-strain response of 8H satin C/Cu- 
alloy composites, with three experimental tensile stress-strain curves for composites with each 
matrix alloy type. Again, the model parameters discussed earlier have been employed. Since 
actual compressive tests were performed using the compressive gripping apparatus, the constraint 
effect is included for the compressive simulations. Clearly, when the same parameters are used 
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Figure 4.11. Comparison of model predictions (V f = 40%) and experimental data for the 
monotonic compressive stress-strain response of 8H satin C/Cu-alloy composites. 
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to model compression as were used to model tension, the actual compressive response of the 
composite is underpredicted. As was the case in tension, the predicted compressive stress-strain 
curves are nearly identical in the elastic region. However, unlike the tensile predictions, the 
compressive predictions are considerably more compliant initially than the experimental stress- 
strain response. This is not because the compressive predictions are more compliant than the 
tensile predictions. In fact, as was shown in Figure 4.9, the initial slope is greater in compression 
than tension due to the constraint effect (in the case of R, = 6x1 O’ 7 , SFT = 70 °F, the initial 
slopes are 11.98 in tension, and 12.10 in compression). Rather, the discrepancy present in 
compression is caused by the noticeably higher initial slopes in the experimental compressive 
stress-strain curves compared to the experimental tensile stress-strain curves. Thus it appears 
that, as modeled, the grip constraint effect is insufficient to account for the differences observed 
in the initial stress-strain behavior of the composite between tension and compression. It is 
likely that a more realistic imperfect fiber-matrix bonding model, which accounts for the closing 
of the interface in compression, would account for some of the discrepancy. Recall that for these 
simulations, the interface is treated as perfectly bonded in the normal direction, but imperfectly 
bonded tangentially. Due to averaging, the tangential debonding does effect the predicted 
normal response of the composite, but it affects the normal tensile and compressive identically. 

Although the envelope of the predicted compressive stress-strain curves falls almost 
completely below the envelope of the experimental compressive stress-strain curves in Figure 
4.11, the correct trends in yielding and hardening are predicted by the model. That is, the 
predicted compressive stress-strain curve for which R, = 6x10’ , SFT = 70 °F, which models the 
C/Cu composite, yields first (and hardens the least), followed by the prediction for R , = 4x1 O’ 7 , 
SFT =120 °F, which models the C/Cu-Ti composite, and finally by the prediction for R, = 
2x1 0’ 7 , SFT =170 °F, which models the C/Cu-Cr composite. Recall that these trends (in the 
experimental stress-strain curves) are opposite to those observed in tension. Since both higher 
values of R , and lower values of stress-free temperature give rise to earlier yielding and lower 
post-yield slope (see Sections 3.3 and 3.4), this trend is expected. Thus, while the quantitative 
correlation between model and experiment in compression is clearly inferior to that in tension, 
the qualitative agreement is superior. 
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Figure 4.12 compares the predicted shear stress-strain response of 8H satin C/Cu-alloy 
composites, with three experimental tensile stress-strain curves for composites with each matrix 
alloy type. Clearly, as was the case in tension, the overall agreement between model and 
experiment is acceptable. Unlike the predicted tensile and compressive stress-strain curves, in 
the elastic region, noticeable differences exist between the three predicted shear stress-strain 
curves. This is because the R, parameter has a much greater effect on the predicted shear 
response of the composite than it does on the predicted tensile response and compressive 
response. In addition, the correct trends in initial modulus, yield stress, and post-yield slope are 
predicted by the model. That is, the predicted shear stress-strain curve for R, = 6x1 O' 7 , SFT = 70 
°F, which models the C/Cu composite, is most compliant overall (i.e. appears lowest in the plot), 
followed by the prediction for R , = 4xl0* 7 , SFT =120 °F, which models the C/Cu-Ti composite, 
and finally by the prediction for R t = 2x1 O' 7 , SFT = 170 °F, which models the C/Cu-Cr 
composite. Recall from Sections 3.3 and 3.4 that while higher R, values give rise to a more 
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Figure 4. 12. Comparison of model predictions (V f = 40%) and experimental data for the 
monotonic shear stress-strain response of 8H satin C/Cu-alloy composites. 
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compliant overall shear response, higher stress-free temperatures cause yielding at lower stresses 
and raise the post-yield slope in shear. Thus, it appears that, for shear, the direct effect of the 
imperfect bonding (modeled via R , ) on yield stress dominates over the effect of the stress-free 
temperature. Since this debonding model provides a level of quantitative as well as qualitative 
correlation in shear that surpasses the correlation in both tension and compression, it can be 
concluded that the imperfect bonding model is reasonably realistic for shear, while it somewhat 
unrealistic for normal loading. It should be noted that, as the strain becomes large, and the 
experimental shear stress-strain curves begin to flatten, the model tends to overpredict the post- 
yield slope of the curves. This may be due to finite local strains or movement of the 
reinforcement weave which cannot be simulated by the model. 

4.3 Cyclic Correlation 

The cyclic tensile, compressive, and shear stress-strain response of 8H satin C/Cu has 
been simulated, as has the combined tension-compression and combined compression-tension 
response. The cyclic response of only one matrix alloy type, Cu-Cr, has been modeled since the 
important effects of matrix alloy type were addressed via simulation of the monotonic response 
of the composite. In order to simulate the 8H satin C/Cu-Cr composite, as before, the parameters 
R , = 2x1 O' 7 and SFT = 170 °F were employed. For cyclic compression, as well as combined 
tension-compression and combined compression-tension, the grip constraint effect was included. 
Other model parameters are identical to those used to generate the monotonic model predictions. 

In the cyclic tension, compression, and shear simulations, unloading for each cycle was 
performed at the identical strain at which unloading occurred in the actual tests being simulated. 
Once the global stress during a given unloading cycle reached zero, the next loading cycle was 
imposed. 

Figure 4.13 compares the predicted cyclic tensile stress-strain response for 8H satin 
C/Cu-Cr with experimental cyclic tension data. The model correlates reasonably well with 
experiment, especially at the lower strains. However, the unloading portions of the stress-strain 
curves indicate clearly that the model overpredicts the modulus of the composite, and the same is 
true for the post-yield slope. Note (see Figure 4.10) that these discrepancies carry over from the 
monotonic correlation. For both monotonic and cyclic loading of C/Cu-Cr, the fact that both the 
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Figure 4.13. Comparison of model predictions (V f = 40%) and experimental data for the cyclic 
tensile stress-strain response of 8H satin C/Cu-Cr composites. 


predicted and experimental stress-strain curves exhibit signs of yielding quite early allows good 
initial correlation. The unloading cycles show the discrepancy in the elastic behavior much more 
clearly. 

At higher strains, the discrepancy between model and experiment is somewhat larger in 
Figure 4.13 than was the case for monotonic tensile loading (Figure 4.10). This is due to scatter 
in the experimental results rather than a major difference in (the envelope of) the model 
predictions. 

The size of the hysteresis loops present in the experimental tensile stress-strain curve is 
underpredicted by the model, indicating that the plastic strain distribution predicted by the model 
is not completely accurate. Strain concentrations that occur in the yam cross-over regions of the 
composite, which are underpredicted by the model, may account for this discrepancy. 

Figure 4.14 compares the predicted cyclic compressive stress-strain response for 8H satin 
C/Cu-Cr with experimental cyclic compression data. The apparent superior model-experiment 
correlation shown in this figure compared to that in Figure 4.1 1 is not due to a difference in the 
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envelope of the model predictions between the two figures. Rather it is due to scatter in the 
experimental results. Though the predicted initial response of the composite is somewhat 
obscured in Figure 4.14, it is in fact noticeable more compliant that the experimental initial 
response. This is similar to the initial correlation for the monotonic compressive response. At 
moderate strains, the envelopes of the predicted and experimental curves agree well, mainly due 
to the fact that the envelope of the experimental curve falls somewhat below the average location 
for C/Cu-Cr compressive stress-strain curves. As was the case in tension, the hysteretic behavior 
of the composite is not well captured by the model, once again indicating a discrepancy in the 
predicted plastic strain distribution for the composite. 

Figure 4.15 compares the predicted cyclic shear stress-strain response for 8H satin C/Cu- 
Cr with experimental cyclic shear data. Agreement between model and experiment is initially 
quite good, but as yielding initiates, the model underpredicts the experimental response of the 
composite somewhat. However, the post-yield shear behavior of the composite is overpredicted, 
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Figure 4.14. Comparison of model predictions (V f = 40%) and experimental data for the cyclic 
compressive stress-strain response of 8H satin C/Cu-Cr composites. 
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Figure 4.15. Comparison of model predictions (V f = 40%) and experimental data for the cyclic 
shear stress-strain response of 8H satin C/Cu-Cr composites. 

and the predicted and experimental stress-strain curves cross. The agreement between the 
(envelopes of the) two curves is thus similar to that shown in Figure 4.12 for monotonic shear 
response of C/Cu-Cr. Note that the scale of the strain axis in Figure 4.15 ends at 0.6 %, while in 
Figure 4.12, it ends at 0.3 %. The hysteresis loops predicted by the model are somewhat larger 
than those predicted for cyclic normal loading, but they are still noticeably smaller than the 
observed hysteresis loops. From an overall quantitative standpoint however, the agreement 
between the predicted and observed cyclic shear response of 8H satin C/Cu is reasonable. 

Figures 4.16 and 4.17 compare model predictions for the combined tension-compression 
and combined compression-tension response of 8H satin C/Cu-Cr with experimental combined 
tension-compression and combined compression-tension results. Note that since the 
experimental combined tests were performed using the compressive gripping apparatus, the 
constraint effect was included in the model simulations. The overall agreement between model 
and experiment in both figures is reasonably good. The envelopes of the predicted and 
experimental response are, for the most part, coincident, and as is the case for the experimental 
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Figure 4.16. Comparison of model predictions (V f = 40%) and experimental data for the 
combined tension-compression stress-strain response of 8H satin C/Cu-Cr composites. 
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Figure 4. 1 7. Comparison of model predictions (V f = 40%) and experimental data for the 
combined compression-tension stress-strain response of 8H satin C/Cu-Cr composites. 
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curves, the predicted curves exhibit a large amount of permanent tensile strain upon unloading 
from tension. Both model and experiment show, however, that upon unloading from 
compression, little permanent compressive strain remains. Since this effect is accurately 
modeled, and the only difference between tension and compression in the model arises from the 
presence of residual stresses, it can be concluded that this effect is caused by the presence of 
residual stresses. 

One difference between the experimental and predicted combined response of the 
composite is that, in the experimental response, the curve for each cycle passes through the point 
of maximum tensile stress reached by the previous cycle. This is not the case for the predicted 
response. In addition, the post-yield slope of the composite for both simulations in both tension 
and compression appear to be somewhat overpredicted. These discrepancies serve as further 
evidence that the plastic strain distributions in the composite predicted by the model are not 
completely accurate. 



5. Summary and Conclusions 


The objective of this investigation was to develop a realistic, versatile, and 
computationally efficient micromechanics model for woven metal matrix composites. The 
model is based on an embedded approach: the original method of cells micromechanics model is 
embedded in the GMC-3D micromechanics model. The nature of the embedded approach, and 
the two micromechanics models used, provided WCGMC with the required versatility. By 
reformulating the equations of GMC-3D, the required computational efficiency was achieved for 
sufficiently refined woven composite microstructures. By including the micro-level features of 
matrix plasticity, porosity, and imperfect fiber-matrix bonding, the model was made realistic for 
a wide range of woven and braided composites. Hence, for the first time woven metal matrix 
composites, with inelastic matrix constitutive behavior, have been modeled herein. 

The model was employed to perform parametric studies of several effects on the 
mechanical response of 8H satin C/Cu. Finally, the model results were compared with 
experimental data for the novel 8H satin C/Cu composite system reported previously by 
Bednarcyk et al. (1997). The major conclusions of this work are summarized in point form 
below: 

• Three distinct geometric models were developed representing the true repeating unit cell for an 
8H satin woven composite. The effect of refining the geometry from its simplest true 
representation to the next level results in an increase in the number of subcells from 192 (TMM 
geometry) to 1280 (CMM and DCMM geometries). Increasing the level of geometric refinement 
gave rise to a more compliant overall predicted tensile and shear response for 8H satin C/Cu, 
with the effect being greatest for the post-yield behavior in shear. 

• Increasing the overall fiber volume fraction of the composite while keeping the yam fiber 
volume fraction constant decreased the predicted tensile, compressive, and shear moduli. This is 
due to the lower transverse tensile modulus and axial shear modulus of the carbon fiber 
compared to the copper matrix which allows the pure copper layers in the unit cell to have a 
greater effect. The predicted tensile, compressive, and shear post-yield behavior, however. 
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became stiffer as the overall fiber volume fraction was increased. This is because the copper 
matrix becomes much more compliant once yielding occurs. 

• Porosity, which is known to exist within the infiltrated fiber yams of the 8-hamess satin 
carbon fiber weave, has a significant effect on the predicted tensile, compressive, and shear 
response of the composite. Increasing the porosity within the infiltrated fiber yams from 0 % to 
7 % to 14 % produced more compliant predicted tensile and shear responses for the composite. 
The simulated effect of porosity was more significant for shear than tension. 

• While the incremental plasticity constitutive theory is realistic for copper at room 
temperatures, due to time-dependent behavior, it is unrealistic for copper at elevated 
temperatures. Through the use of a previously developed version of lamination theory utilizing a 
viscoplastic constitutive model, it was determined that to realistically estimate the residual 
stresses in 8H satin C/Cu using incremental plasticity, an effective stress-free temperature (SFT) 
below 250 °F should be used in WCGMC. 

• For SFT values below 250 °F, increasing the SFT caused yielding at lower stresses and higher 
post-yield slopes during subsequent simulated mechanical loading. The predicted yield stresses 
were lower because subcells with higher residual stress tend to yield at lower applied global 
strains. The increased post-yield slopes for higher SFTs were caused by the higher degree of 
non-uniformity in the residual stress state in the composite. This non-uniformity implies that 
many subcells have stress states which are not favorable to the accumulation of additional plastic 
strain. In fact, many subcells which yield during the simulated cool-down load elastically upon 
subsequent simulated mechanical loading. Differences in the magnitude of the SFT effect are 
apparent because, in most subcells, the residual stress state is tensile. Thus, in the case of 
simulated tension, more subcells load plastically from the start. 

• Imperfect fiber-matrix bonding has a significant effect on the predicted tensile and shear 
response of the composite. While the predicted in-plane shear response of the composite was 
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affected by onJy the tangential debonding parameter, the predicted in-plane tensile response was 
affected by both the normal and tangential debonding parameters. The effect of the tangential 
debonding parameter on the predicted normal response of 8H satin C/Cu is caused mainly by the 
averaging procedure employed in the original method of cells. 

• Realistic model parameters were chosen in order to obtain the best possible agreement 
between model and experiment. Due to the unrealistic nature of the normal imperfect bonding 
model in the presence of interfacial compression, the normal debonding parameter was not used. 
Different values of the tangential debonding parameter were chosen to simulate each matrix alloy 
type, and to provide good agreement of the model predictions with experiment for shear. To 
simulate the varying effect of longitudinal stress-dependent sliding at the interface during cool- 
down for each matrix alloy type, different SFT values were employed for each matrix alloy type 
as well. 

• Utilizing the above parameters in WCGMC provided good quantitative correlation between 
model and experiment for monotonic tension and shear. The model predictions for monotonic 
compression were somewhat more compliant than the experimental compressive results which 
indicates the need for a superior model for fiber-matrix debonding. The qualitative trends based 
on matrix alloy type (which correspond to trends in fiber matrix bonding) were correctly 
predicted for monotonic compression and shear, and, most notably, for monotonic tension as 
well. In tension, while the trend caused by the chosen values of the tangential debonding 
parameter was opposite to the trend observed in experiment, the trend caused by the chosen SFTs 
is the same as the observed trend. Since the effect of SFT choice dominates the predicted tensile 
response, the model correctly predicts the trend observed in experiment. 

• The predictions of WCGMC correlate reasonably well with experiment for cyclic tension, 
compression, and shear, as well as for combined tension-compression and combined 
compression-tension. The main discrepancy between the model predictions and experiment in 
these cases involves the underprediction of the observed hysteresis loops. This discrepancy 
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indicates that the plastic strain fields predicted by the model may be somewhat inaccurate and 
that a more realistic debonding model is needed. 
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where Rp = {$ + N 2 + N a 
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x i r) =/Ya >, 

1 (VA^A) 


N, 

0,2 = 2d/Y G, , , 

/=I 

N J 

tY = 2dhTG, , 

12 ^ (V + *2+^+^) 

^=-| G <vA 




where R y = y + N 2 + N a + , 


where G^j are the components of the inverse of the global matrix G in equation (4.37), f* are the components 
of the thermal vector f 4 , and f p are the components of the plastic vector f p . 
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